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Schematic  diagram  of  the  experimental  procedure  for 
obtaining  bulk  measurements 


Schematic  diagram  of  the  experimental  arrangement  for 
obtaining  bulk  measurements  


The  inside  of  the  vacuum  chamber  is  shown 


The  vacuum  chamber  is  shown  on  theleft  and  the  inter- 
ferometer with  collection  optics  and  detector  on  the  right; 
the  beam  displacer  is  mounted  in  between  for  easy 
alignment  


Com 

fo’- 


reflectance  as  a function  of  incidence  angle 
different  values  of  complex  refractive  index 


<_uvity  spectrum  (solid  curve)  of  carbon 
ite)  from  400  cm”^  (25  pm)  to  2000  cm- 1 (5.0  pm) 
measure  at  room  temperature  and  an  angle  of  incidence 
of  45  degrees 

The  index  of  refraction  (n)  of  carbon,  from  400  cm"  ^ 
to  2000  cm"  as  determined  from  a Kramers-Kronig 
analysis  of  the  experimental  reflectivity  data  is  shown 
as  the  solid  curve 


The  extinction  coefficient  (k)  of  carbon  from  400  cm" 
to  2000  cm-1  as  determined  from  a Kramers-Kronig 
analysis  of  the  experimental  reflectivity  data  is  shown 
as  the  solid  curve 


The  index  of  refraction  of  carbon  materials  in  the 
wavenumber  region  400  cm"  ^ to  2200  cm- 1 . . . . 


The  extinction  coefficient  of  carbon  materials  in  the 
wavenumber  region  400  cm- ^ to  2200  cm”  ^ 


The  reflectivity  spectrum  of  MgO  in  the  400  cm"^  to 
2000  cm"^  region  


The  reflectivity  spectrum  of  MgO  in  the  400  cm"*  to 
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SECTION  I 


INTRODUCTION 


Scattering  and  thermal  emission  of  electromagnetic  radiation 
from  particles  of  various  compositions  have  taken  on  significance  in 
several  areas  of  technology  during  the  past  few  years.  In  addition  to 
the  importance  of  these  processes  in  certain  environmental  problems, 
emission^  absorption  and  scattering  of  radiation  from  particles  plays  a 
crucial  role  in  certain  Air  Force  Systems  applications.  The  inter- 
actions of  particles  with  energy  sources  can  significantly  change  the 
electromagnetic  radiation  field  from  that  expected  in  the  absence  of 
particles. 

To  our  knowledge,  at  the  beginning  of  the  present  study  there 
were  no  measurements  of  single  particle  emissivity  in  the  infrared  at 
room  temperature  or  elevated  temperatures.  To  compound  the  pro- 
blem of  a lack  of  experimental  data,  there  were  apparent  discrepan- 
cies between  visible  emissivity  measurements  based  on  bulk  and  on 
particulate  matter.  As  a consequence,  the  reliability  of  using  bulk 
refractive  index  data  to  calculate  particle  emissivity  in  the  infrared 
was  open  to  question.  Furthermore,  there  was  a dearth  of  high  tem- 
perature bulk  refractive  index  data  for  materials  of  interest. 

The  Particle  Optical  Properties  Measurements  (POPM)  Program 
was  proposed  to,  and  jointly  funded  by  AFRPL  and  SAMSO/ABRES,  to 
attempt  to  provide  answers  to  these  questions.  The  program  agreed 
upon  between  AFRPL,  SAMSO/ABRES,  and  The  Aerospace  Corporation 
in  June  of  1974  consisted  of  four  tasks: 

(1)  To  perform  an  exhaustive  and  critical  compilation  of  exist- 
ing measurements  of  bulk  refractive  index  at  room  temp- 
erature and  above  for  four  materials:  Carbon,  MgO,  ZrO^ 
and  Al^O^. 

To  fabricate  an  apparatus  and  procedures  for  the  deter- 
mination of  the  infrared  complex  refractive  index  of  "bulk" 


(2) 


samples  of  materials  at  high  temperatures.  (By  ''bulk'' 
samples  we  mean  samples  which  have  dimensions  such 
that  a polished  surface  ~ 3.  2 X 7.  6cm  can  be  prepared  on 
the  sample.  ) This  apparatus  was  then  to  be  used  to 
measure  high  temperature  infrared  complex  refractive 
indices  for  carbon,  MgO,  ZrO^  and  Al^O^. 

(3)  To  design  and  construct  a device  to  measure  the  emissivity 
of  single  spherical  particles  heated  to  high  temperature 

by  laser  radiation.  This  apparatus  was  to  be  used  to 
determine  the  emissivity  of  single  particles  of  various 
sizes  in  the  micron  range  formed  from  the  four  materials: 
carbon,  MgO,  ZrO^  and  Al^O^.  The  observed  emissivi- 
ties  were  to  be  compared  with  the  emissivities  computed 
from  the  high  temperature  bulk  refractive  indices. 

(4)  To  use  the  apparatus  built  for  Task  3 to  measure  the 
emissivity  of  particles  provided  by  AFRPL. 

The  first  of  these  tasks  has  been  completed  and  the  results  pub- 
lished in  a combined  AFRPL/Aerospace  report.  ^ ^ This  task 

will  not  be  discussed  further  in  this  report.  The  fourth  task  was 
deleted  by  mutual  agreement  when  sufficient  quantities  of  particles  to 
make  both  bulk  and  particle  measurements  were  not  available. 

Tasks  2 and  particularly  3 involve  combinations  of  state-of-the- 
art  measurement  techniques.  This  report  discusses  the  procedures 
developed  and  results  obtained  on  these  two  tasks.  The  second  section 
discusses  the  bulk  refractive  index  experiment  and  results.  The  third 
section  provides  the  theory  and  computed  particle  emissivity  based  on 
Mie  theory  and  the  refractive  index  results  presented  in  Section  II. 
Section  IV  discusses  the  techniques  and  apparatus  developed  for  the 
single  particle  emissivity  experiment. 


SECTION  II 


HIGH  TEMPERATURE  BULK  REFRACTIVE  INDEX  EXPERIMENT 


A.  EXPERIMENTAL  DESCRIPTION 

The  optical  constants,  n and  k,  of  the  materials  of  interest  were 
obtained  in  the  wavelength  range  4.  5ijjm  to  25.0^i,m,  and  the  measure- 

O O 

ments  were  made  with  sample  temperatures  of  300  K and  ~ 600  K for 
all  the  materials  except  carbon,  which  was  measured  at  300 °K  only. 

The  complex  index  of  refraction  was  extracted  from  the  reflectance 
spectrum  by  applying  what  has  come  to  be  known  as  Kramers-Kronig 
analysis  and,  when  necessary,  the  classical  dispersion  model  calcula- 
tion. 

In  1952  Robinson^^^^'  first  applied  the  K-K  analysis  to 

reflectance  spectra  that  were  obtained  from  normally  incident  radia- 
tion. Since  then  the  technique  has  been  extended  to  make  it  applicable 
for  all  angles  of  incidence  between  0 and  90  degrees,  ^ 

provided  linearly  polarized  light  is  used.  The  classical  dispersion 
model  is  also  a time-honored  technique^^^^^*  and  can  be 

applied  to  reflectance  spectra  obtained  with  any  polarization  and  angle 
c>f  incidence.  These  two  theories  are  presented  in  detail  in  Section  Il-B 
of  this  report. 

The  experimental  procedure,  shown  schematically  in  Figure  II- 1, 
for  measuring  the  reflectance  spectrum  is  to  chop  the  radiation  flux  of 
an  infrared  source,  collimate  it,  linearly  polarize  the  light  with  the 
electric  field  perpendicular  to  the  plane  of  incidence,  reflect  the  radia- 
tion from  the  surface  of  the  sample,  spectrally  modulate  the  light  with  a 
fast  scanning  Michelson  interferometer,  and  finally,  collect  the  radiation. 

In  this  investigation  ~ 1600  interferograms  were  obtained  and 
coherently  added  to  improve  the  signal-to-noise  ratio.  For  each 
sample^ calibration  spectra  were  obtained.  The  calibration  standard 
was  a gold  mirror.  The  sequence  of  measurement  was:  Run  calibra- 
tion spectra;  run  sample  spectra;  run  calibration  spectra. 


VACUUM  TANK 


Figure  II- 1.  Schematic  diagram  of  the  experimental  procedure  for  obtaining  bulk  measurements. 
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Experimental  Layout 


The  experimental  layout  for  the  bulk  Index  of  refraction  meas- 
urements Is  shown  In  Figure  II-Z.  The  path  of  radiation  within  the 
tank  Is  shown  In  Figure  II-3.  The  Infrared  radiation  Is  provided  by 
an  Opperman  source,  Perkln-Elmer  457-0244,  which  Is  a ceramic 
coated  resistance  wire.  The  radiation  Is  chopped  by  a small  mirror, 

0.  2 X 0.3cm,  mounted  on  a torsional  tuning  fork,  Bulova  type  L45, 
tuned  to  lOkHz,  having  an  angular  amplitude  of  two  degrees  peak  to 
peak  and  which  has  been  placed  close  to  the  source  exit  aperture. 

The  light  Is  then  rendered  parallel  by  a 31cm  focal  length  spherical 
mirror  used  slightly  off-axis.  The  parallel  beam  from  the  source  Is 
tj  en  folded  and  passed  through  a Perkln-Elmer  wire  grid  polarizer, 
a gold  wire  grid  deposited  on  silver  bromide.  The  efficiency  of  this 
polarizer  Is  about  98  percent  at  2pm  and  Increases  to  more  than  99.  8 
percent  at  lOpm;  the  transmission  Is  20  percent. 

The  source  beam,  polarized  with  Its  electric  field  perpendicular 
to  the  plane  of  Incidence,  and  with  a diameter  of  2.  2cm,  strikes  the 
surface  of  the  sample  at  an  angle  of  Zrn  relative  to  the  optical  axis  of 
the  Interferometer  and  Is  reflected  along  that  axis;  the  angle  of  Inci- 
dence Is  O'.  The  equipment  Is  designed  to  provide  a range  In  a 
between  25  and  80  degrees  with  some  vignetting  at  the  very  largest 
angles.  The  angle  a Is  measured  with  a digital  angular  encoder 
(Norden  Model  5-36  BCDE360L)  to  an  accuracy  of  ~ 0.  05  degrees. 

For  the  measurements  presented  In  this  report  a was  set  at  45  degrees. 

The  sample  dimensions  are  7.  6 X 3.  2 x 1.  3cm.  The  front  and 
rear  surfaces  are  cut  so  that  they  are  not  parallel.  This  Insures  that 
the  radiation  reflected  from  the  back  surface  does  not  enter  the  Inter- 
ferometer. The  sample  has  an  electrical  heater  pressed  against  the 
rear  surface.  The  mount  for  the  sample  and  heater  Is  designed  to  serve 
as  a heat  shield;  It  Is  made  of  tantalum  and  surrounds  all  but  the  front 
of  the  sample.  To  minimize  heat  losses  due  to  conduction,  the  sample 

Is  mounted  with  low-conductivity  ceran^lc  parts;  convective 
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Figure  11-3.  The  inside  of  the  vacuum  chairiber  is  shown.  A laser  beam  traces  tlie  optical  path 
from  the  source,  lower  right  of  the  photograph,  to  the  exit  froni  the  chamber, 
upper  left.  The  flexible  tubing  is  used  to  fill  and  vent  the  liquid  nit»'Ogen  reservoir 
at  the  bottoni  of  the  photograph,  which  cools  Uie  source  optics. 
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heat  losses  are  eliminated  by  evacuating  the  whole  system.  Surface 

O 

temperatures  of  up  to  700  K have  been  reached  by  this  method;  attempts 
for  higher  temperatures  resulted  in  heater  failure.  A new  heater  sys- 
tem is  now  being  tested  which  is  designed  to  allow  the  attainment  of 

O 

surface  temperatures  up  to  and  in  excess  of  1000  K. 

An  array  of  spring-loaded  thermocouples  touching  the  front  sur- 
face of  the  sample  was  used  to  detect  thermal  gradients.  The  largest 
difference  observed  on  the  surface  was  about  four  degrees  at  600 °K. 
During  data  acquisition  the  thermocouple  array  is  not  used;  the  tem- 
perature is  then  monitored  with  a thermocouple  touching  the  surface 
out  of  the  field-of- view  of  the  interferometer. 

After  reflection  from  the  sample,  the  source  beam  leaves  the 
vacuum  chamber  through  an  Irtran  6 window.  Light  emerging  from 
the  vacuum  chamber  is  brought  into  the  interferometer  by  means  of 
a beam  displacer  inserted  in  the  optical  train  for  convenience  in 
alignment.  These  optical  components  outside  the  vacuum  tank  are 
shown  in  the  left  of  Figure  II-4. 

The  Digilab  296  fast  scan  Fourier  transform  interferometer, 
shown  in  the  center  of  Figure  II-4,  is  of  the  Michelson  type  modified 
so  that  double-sided  interfe rograms  can  be  obtained.  Since  double- 
sided interferograms  do  not  require  phase  correction,  the  possibility 
of  errors  being  introduced  by  a phase  correction  is  eliminated.  The 
radiation  leaving  the  interferometer  is  collected  by  a 10cm  diameter, 
76cm  focal  length  mirror  and  is  detected  by  a liquid  helium-cooled, 

copper  doped  germanium  detector,  Raytheon  QKN1546.  The  detector 

-3  2 

characteristics  are:  KRS-5  entrance  window,  7.85  X 10  cm  area, 
five  degrees  field-of-view,  6.  8 X 10  N.  E.  P.  To  limit  the 

bandpass,  a cold  filter.  Spectrum  Systems  54958X,  was  mounted  inside 
the  dewar.  The  transmission  of  the  filter  is  greater  than  60  percent 
in  the  2.8ym  to  22.  0y,m  range  and  less  than  0.5  percent  from  0.4|j,m 
to  2.  0|im. 
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The  source  housing,  collimating  and  folding  optics  are  liquid 
nitrogen-cooled  to  lower  the  unwanted  thermal  background  radiation, 
which  might  overload  the  detector.  In  addition,  the  chamber  walls 
are  shielded  by  a liquid  nitrogen-cooled  shroud.  Near  the  end  of  the 
experimental  measurements  the  vacuum  integrity  of  this  chamber 
shroud  failed.  However,  little  degradation  of  experimental  perform- 
ance was  observed  after  this  failure. 

2.  Signal  Processing 

The  interferogram  was  sampled  when  the  optical  path  difference 
between  the  stationary  mirror  and  moving  mirror  was  increased  by 

_4 

0.63299  X 10  cm,  the  wavelength  of  the  He-Ne  laser  line;  for  this 
rate  of  data  acquisition  sampling  theory  dictates  that  the  largest  wave- 
number  with  non-zero  intensity  should  be  ~ 7900cm~^. 

The  scan  velocity  of  the  interferometer  mirror  is  ~0.5cm/sec, 
therefore,  our  largest  wavenumber,  7900cm~\  results  in  a frequency 
of  3950Hz.  (This  means  that  a monochromatic  beam  with  v = 7900cm~ 
would  have  as  its  interferogram  a cosine  wave  of  frequency  3950Hz.  ) 

The  signal  that  reaches  the  detector  is  a lOkHz  carrier  wave 
(due  to  the  chopper)  modulated  by  lower  frequencies  (due  to  the  inter- 
ferometer) with  the  highest  of  these  being  ~ 3950Hz.  This  signal  is 
amplified  and  then  synchronously  detected  with  a Princeton  Applied 
Research  lock-in  amplifier,  model  124A. 

The  last  stage  of  the  lock-in  amplifier,  which  provides  the  final 
filtering  of  the  demodulated  signal,  has  a high  frequency  cutoff  less 
than  3950Hz  and  so  it  is  not  used.  Instead,  a bandpass  filter  with  a 
6kHz  cutoff  frequency  and  an  attenuation  of  -40dB  at  lOkHz  has  been 
built  from  commercial  active  bandpass  filter  components  and  is  used 
to  filter  the  demodulated  signal  before  being  fed  to  the  data  handling 
system. 
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The  data  handling  system  samples  the  filtered  analog  signal 

from  the  lock-in  amplifier  when  commanded  by  the  interferometer. 

The  sampled  signal  is  converted  to  a 12-bit  digital  number  in  an 

analog  to  digital  converter.  The  resulting  number  is  recorded  on 

computer  compatible  tape  for  subsequent  Fourier  transformation  and 

processing  of  the  spectra  by  the  Aerospace  Control  Data  7600  compu- 

r II  - 9l 

ters  using  standard  techniques.*^  * ' ^ 

3.  Samples 

Bulk  samples  had  to  be  at  least  3.  2 x 7.  6 x 1. 3cm,  of  high 
purity  and,  when  appropriate,  single  crystals.  Since  this  last  condi- 
tion could  not  be  met  for  all  of  the  materials,  sources  and  methods 
were  sought  that  could  produce  bulk  samples  from  pure  powders.  The 
surfaces  of  the  samples  were  polished  with  a compound  made  of  0.  3u,m 
AI2O2  particles  and  distilled  water,  and  Linde  type  A polishing  com- 
pound. They  were  then  cleaned  with  an  ultrasonic  cleaner  and  different 
solvents;  acetone  was  first  used,  then  a detergent  and  ammonia  com- 
bination, and  finally,  distilled  water.  Impurity  contamination  of  the 
samples  is  felt  to  be  negligible,  however,  the  physical  process  of  fx>l- 
ishing  and  cleaning  can  affect  the  surface  of  the  samples  (i.  e.  , intro- 
duction of  strains)  and  thus  may  affect  the  reflectivity  measurements. 
This  is  discussed  in  more  detail  in  the  section  on  the  experimental 
observations  on  Al^O^. 

a.  Carbon  (C) 

Single  crystals  of  graphite  can  be  produced  synthetically. 
However,  available  sizes  are  not  large  enough  to  use  in  this  experiment. 
The  sample  we  have  used  is  CEP  graphite,  produced  by  Union  Carbide. 
This  is  a mixture  of  lamp  black  and  binder  pressed  and  baked  at  a very 
high  temperature. 
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b.  Magnesia  (MgO) 

A single  crystal  of  MgO  was  obtained  from  the  ROC/RIC 
Corporation  of  Sun  Valley,  California.  This  material  is  99.  99  percent 
pure. 


c.  Zirconia  (ZrO^) 

Single  crystals  of  zirconium  dioxide,  large  enough  for  our 
needs,  do  not  occur  in  nature  and  are  not  synthetically  produced.  The 
only  available  choice  was  a polycrystalline  piece  formed  by  hot  pres- 
sing Zr02  powder. 

The  Cerac  Corporation  supplied  Aerospace  with  a hot  pressed 
ZrO^  sample  that  is  98  percent  ZrO^  plus  two  percent  HfO^,  an 
impurity  which  cannot  be  economically  removed  from  Zr02. 

d.  Alumina  (Al^O^) 

Synthetic  pure  crystals  of  alumina  were  found  to  be  avail- 
able. The  Union  Carbide  Corporation  supplied  us  with  a single  crystal 
of  pure  Al^O^* 


B.  DATA  REDUCTION  METHODS 

1,  Angular  Reflectance  Analysis,  Avery  Method 

The  index  of  refraction  and  the  extinction  coefficient  can  be 
determined  from  a measurement  of  the  reflectance  Rj^^  and/or  R^^ 
(radiant  flux  polarized  parallel  and  perpendicular  to  the  plane  of  inci- 
dence, respectively)  by  the  method  of  Avery, 

This  method  is  based  on  the  generalized  Fresnel  reflectance 

V u u [Ref.  II-ll] 

equations  which  can  be  written*^ 
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R 


D - (Q  ~ cos  6)  + P 

^1  “ 2 2 
(Q  + cos  0)  + P 


and 


11 


= R 


(Q  - sin  0 tan  0)^  + ^ 

■*■  (Q  + sin  0 tan  0)  + P 


- P^  = n^  - - sin^  0,  and 


QP  = nk. 


(II-l) 


(II-2) 


k- 

r 


where  n is  the  index  of  refraction  (real  part),  k is  the  extinction  coef- 
ficient, and  0 is  the  angle  of  incidence.  Thus  in  principle,  a measure- 
ment of  Rj^  and/or  R^^  as  a function  of  0 will  allow  determination  of 

n and  k.  A direct  method  for  inverting  these  equations  to  determine 

j • .u  1-.  * [Refs.  II-ll  -11-13] 

n and  k has  been  discussed  recently  in  the  literature. 

Unfortunately,  when  the  reflectance  data  for  Al^O^  were  analyzed  by 

this  method,  the  equations  proved  to  be  ill-conditioned  under  certain 

conditions  and  unreasonable  values  of  n and/or  k were  sometimes 

obtained.  Invariably  in  these  situations,  small  changes  in  the  input 

reflectance  data  produced  large  changes  in  n and/or  k.  The  reason  for 

this  ill-conditioning  is  not  too  difficult  to  discover.  For  instance,  if 

n is  small  1)  and  k large  7),  both  Rj^  and  R^^  are  > 0.8,  and  change 

little  with  0 or  increasing  k.  Another  example  of  ill-conditioning  is 

illustrated  in  Figure  II-5.  Plotted  is  Rj^  as  a function  of  the  angle  of 

incidence  0 = 45  degrees  for  the  three  cases:  n = 3.  5,  k = 3.  0;  n = 4.  5, 

k = 3.  0;  and  n = 5.  5,  k = 2.  0.  The  curves  do  not  appreciably  differ 

from  each  other  until  0 is  greater  than  about  72  degrees.  Beyond  0 = 75 

degrees  the  curves  are  significantly  different,  however,  measurements 

at  these  large  angles  of  incidence  require  very  large  samples  (because 

of  the  foreshortening  effect),  and  very  accurate  angular  measurements. 
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due  to  the  rapid  change  of  Rj ^ with  respect  to  0 for  6 ^ 75  degrees. 
Since  neither  of  these  conditions  are  readily  achieved  in  our  experi- 
mental apparatus,  we  have  employed  either  a Kramer s-Kronig  analy- 
sis or  a classical  dispersion  analysis  to  obtain  the  optical  constants. 


2. 


Kramers-Kronig  Analysis 

The  amplitude,  Rj^,  and  phase,  of  the  complex  reflectivity 
coefficient  are  not  independent  but  can,  in  general,  be  related  by  a 
Kramers-Kronig  dispersion  relation.* 


V + v„ 
c 

d 

, Ir,  / vl  1/2 

V - V 

c 

dv 

ln|Rj^(v)| 

dv 


(II- 3) 


where  v is  the  wavenumber  of  the  reflected  wave.  The  generalized 
Fresnel  equations  can  then  be  written  in  terms  of  this  amplitude  and 
phase. 


(1  - R,  ) cos  0 


Q = 


P = 


1 + Rj^  - 2 \/  R^  cos  0 
-a/i:  cos  0 sin  ij) 


(II -4) 


1 + R^  - 2 /r^  COS  if) 


where  the  quantities  have  their  previous  definitions.**  Finally,  n and 
k are  obtained  from  P and  Q by  solving  equation  11-2. 


i 


i ; 


*Kramers-Kronig  analyses  have  been  treated  extensively  in  the  liter- 
ature. The  analysis  outlined  here  follows  Ref.  11-2. 

**Similar  but  more  complicated  equations  can  be  written  for  R^^* 
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The  disadvantage  of  this  analysis  Is  that  0 depends  on  the  behav- 
ior of  Rj^(v)  over  all  of  wavenumber  space.  In  practice,  the  reflectiv- 
ity Is  measured  over  only  a finite  portion  of  wavenumber  space  and 
realistic  approximations  for  the  behavior  of  R^(v)  outside  this  region 
are  required  for  accurate  evaluation  of  the  Integral. 

Most  of  the  methods  discussed  In  the  literature  assume  some 
simple  functional  form  for  the  reflectivity  In  the  outside  region.  In 
some  cases  the  form  of  the  function  Is  justified  by  a theoretical  model 
and  In  other  cases  It  Is  just  a simple  function  which  behaves  In  a 
physically  reasonable  manner.  11-16] 

We  did  not  use  any  special  functional  form  to  extrapolate  the 
reflectivity  curves,  but  rather  tried  to  Incorporate  as  much  as  possi- 
ble any  published^ experimental  data  In  the  outside  regions. 
In  some  cases  the  results  of  classical  dispersion  theory  were  Incor- 
porated and  this  Is  discussed  In  the  next  section.  The  reflectivity 
was  then  extrapolated  graphically  to  fit  these  data  points  and  to  join 
them  smoothly  to  our  measured  reflectivity  spectra. 

This  curve  was  then  digitized  for  computer  calculation  of  the 
phase  by  the  Kramers-Kronlg  Integral.  The  amplitude  and  phase  of 
the  reflectance  were  used  to  obtain  P and  Q (equation  II-4)  and  finally 
n and  k determined  from  solutions  of  equation  II-2. 


3.  Classical  Dispersion  Analysis 

The  classical  dispersion  analysis  of  reflectivity  data  Is  based 
on  a model  which  approximates  the  material  as  a system  of  damped 
harmonic  oscillators.  Each  oscillator  Is  described  by  Its  strength 
Py  width  Yj,  and  frequency  Vj.  The  expression  for  the  complex 
dielectric  constant  (e)  from  tills  model  Is 


e(v)  = ^ 


j 1 - (v/Vj)  - l(v/Vj)  (Yj/Vj) 
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where  is  the  high  frequency  limit.  The  Index  of  refraction,  n,  and 
the  coefficient  of  extinction,  k,  are  related  to  e by 


n + Ik  = y e 


(11-6) 


The  reflectivity  can  then  be  calculated  from  Fresnel's  equation  if  the 
constants  in  equation  II-5  are  known  or  assumed.  These  constants 
are  called  the  dispersion  parameters  and  the  process  of  determining 
them  is  dispersion  analysis.  The  parameters  are  varied  until  agree- 
ment is  obtained  with  experiment.  The  criteria  used  was  to  minimize 
the  sum  of  the  square  errors.  This  was  accomplished  with  a non-linear 
least  squares  fitting  computer  program  which  systematically  varied  the 
parameters  until  a best  fit  was  obtained. 


C.  EXPERIMENTAL  RESULTS 

The  Kramers-Kronig  analysis  was  used  to  reduce  the  data  for 
carbon,  Al^O^,  and  ZrO^  where  we  felt  we  could  adequately  extrapo- 
late the  reflectance  data  to  the  lower  (v  < 400cm"  and  the  higher 
(v>  2200cm  wavenumber  regions.  The  MgO  data  was  reduced 
using  the  classical  dispersion  analysis  method  where  we  felt  we 
could  not  safely  extrapolate  the  reflectance  data  to  the  lower  (v  < 
400cm  wavenumber  region.  In  general,  we  feel  the  Kramers- 
Kronig  method  of  analysis  is  superior  to  the  classical  dispersion 
method  of  analysis,  since  the  former  is  independent  of  any  model 
assumptions  concerning  the  behavior  of  the  optical  parameters 
(n  and  k)  of  the  material,  where  as  the  latter  assumes  the  behavior 
of  the  optical  parameters  of  the  material  may  be  modeled  as  a sys- 
tem of  damped  harmonic  oscillators. 


1 . Carbon 

The  room  temperature  reflectivity  spectrum  of  a microcrystal- 
line,  polished  slab  of  graphite,  at  an  angle  of  incidence  of  45  degrees, 
is  shown  in  Figure  II-6.  The  resolution  is  ~ 15cm" 

The  observed  reflectance  was  smoothly  extrapolated  below  400cm" \ 
as  shown  in  Figure  II-6,  and  a Kramers -Kronig  analysis  was  performed. 
The  resulting  values  of  n and  k are  shown  in  Figures  II-7  and  11-8, 
respectively. 

The  variations  of  the  observed  reflectivity  about  the  generally 
smooth  trend  of  the  graphite  reflectivity  spectra  are  consistent  with 
the  expected  noise  level.  It  was  discovered,  empirically,  that  a 
curve  of  the  form 

= A + Bv*^  (II-7) 

could  very  accurately  fit  our  experimental  data.  This  smoothed 
reflectance  spectrum  is  shown  as  the  dashed  curve  in  Figure  II-6.  A 
Kramers-Kronig  analysis  was  performed  and  the  resulting  n and  k are 
shown  as  the  dashed  curves  in  Figures  II-7  and  II-8,  respectively. 

For  some  situations  the  scaling  in  Figures  II-7  and  II-8  may  be 
too  coarse,  and  therefore  n and  k (as  deduced  from  the  experimental 
data  and  the  smoothed  data)  are  given  as  a function  of  wavenumber,  in 
12cm  ^ increments,  in  Table  II-l,  from  v = 400cm"^  to  2200cm"^. 

In  an  effort  to  ascertain  the  experimental  errors  in  our  meas- 
urements, and  to  estimate  the  resulting  errors  in  n and  k,  the  follow- 
ing analysis  was  performed.  At  ~ 12cm” intervals  the  quantity 
A = I Rj^(s)  - Rj^(e)|  was  calculated  where  Rj^(e)  is  the  experimentally 
measured  reflectance  (Figure  U-6),  and  Rj^(s)  is  the  smoothed  reflec- 
tance from  Figure  II-6.  These  As  were  then  used  to  estimate  the 
error  made  in  measuring  reflectivity  in  two  hundred  wavenumber 
intervals  from  400cm"^  to  2200cm"  The  results  are  given  in  Table 
II-2.  The  ARs,  the  reflectivity  errors,  were  calculated  in  a root- 
mean-square  (rms)  fashion,  that  is 
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Figure  H-7.  The  index  of  refraction  (n)  of  carbon,  from  400cm  to  2000cm  as  determined 

from  a Kramers-Kronig  analysis  of  the  experimental  reflectivity  data  is  shown  as 
the  solid  curve.  The  dashed  curve  is  n as  determined  from  a Kramers-Kronig 
analysis  of  the  smoothed  reflectivity  data. 


Table  II-l.  Measured  Index  of  Refraction  and  Extinction  Coefficient 


for  Carbon  from  v = 400cm 


to  2200cm 


(cm  ) 

*^smoothed 

'^expe  r imental 

smoothed 

experimental 

400 

2.  75 

2.  92 

3.  36 

3.  32 

12 

2.  78 

3.  07 

3.  28 

3.43 

24 

2.  81 

3.  00 

3.  19 

3.  64 

36 

2.83 

2.  50 

3.  11 

3.48 

48 

2.  85 

2.  64 

3.  03 

3.  05 

60 

2.  88 

2.  66 

2.  96 

3.  09 

72 

2.  90 

2.  55 

2.  89 

2.  76 

84 

2.  91 

2.  77 

2.  82 

2.  51 

96 

2.  93 

3.  13 

2.  76 

2.  40 

508 

2.  95 

3.  21 

2.  70 

2.  76 

20 

2.  96 

3.01 

2.  65 

62 

32 

2.  97 

3.  10 

2.  60 

2.  52 

44 

2.  98 

3.  26 

2.  55 

2.  55 

56 

2,  98 

3.06 

2.  50 

2.  76 

68 

2.  99 

2.  99 

2.45 

2.  38 

80 

3.  00 

3.  36 

2.  40 

2.  49 

92 

3.  01 

3.  20 

2.  36 

2.  75 

604 

3.02 

2.  78 

2.  31 

2.  73 

16 

3.  03 

2.  76 

2.  28 

2.  27 

28 

3.  03 

3.  00 

2.  24 

2.  23 

40 

3.  03 

3.08 

2.  21 

2.  30 

52 

3.  04 

2.  94 

2.  17 

2.  35 

64 

3.  04 

2.  86 

2.  13 

2.  23 

76 

3.  04 

2.86 

2.  10 

2.  05 

88 

3.  05 

3.  00 

2.  07 

2.  00 

. I 
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Table  II- 1.  (Continued) 


00 

12 

24 

36 


smoothed 

’^experimental 

smoothed 

^expe  r imental 

3.07 

3.  04 

1.  50 

1.  50 

3.  07 

3.  12 

1.49 

1.44 

3.07 

3.  15 

1.47 

1.45 

3.  07 

3.  18 

1.46 

1.  44 

3.  07 

3.  20 

1.45 

1.48 

3.  07 

3.  16 

1.43 

1. 48 

3.06 

3.  17 

1.42 

1.46 

3.  06 

3.  20 

1.  40 

1.48 

3.  06 

3.  15 

1. 39 

1.  57 

3.  07 

3.  04 

1.  38 

1.48 

3.  06 

3.  10 

1. 36 

1.43  ■ 

3.07 

3.  07 

1.  35 

1.44 

4.  06 

3.  06 

1. 34 

1.  40 

3.  06 

3.  06 

1. 33 

1.  36 

3.  06 

3.09 

1.  32 

1.  35 

3.  06 

3.  09 

1.  30 

1.  35 

3.  06 

3.  09 

1.  29 

1.  35 

3.  06 

3.  08 

1. 28 

1.  33 

3.06 

3.  09 

1.  27 

1.  33 

3.  06 

3.  07 

1. 26 

1.  30 

3.  06 

3.  13 

1.  25 

1.  28 

3.  06 

3.  10 

1.  24 

1.  36 

3.  05 

3.  04 

1.  23 

1.29 

3.  06 

3.09 

1.  22 

1.  26 

3.  06 

3.  12 

1.21 

1.  28 
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Table  II-l . (Continued) 


V 

(cm  ) 

^smoothed 

^experimental 

smoothed 

^experimental 

1300 

3.  06 

3.08 

1.21 

1.  31 

12 

3.  05 

3.  08 

1.  20 

1.  24 

24 

3.  05 

3.  14 

1.  19 

1.  28 

36 

3.  05 

3.  08 

1.  18 

1.  33 

48 

3.  05 

3.  03 

1.  17 

1.  32 

60 

3.  05 

2.  98 

1.  16 

1.  29 

72 

3.  05 

2.  96 

1.  15 

1.  19 

84 

3.  05 

3.07 

1.  14 

1.  18 

96 

3.  05 

3.  05 

1.  13 

1.  24 

1408 

3.  05 

3.  01 

1.  13 

1.  25 

20 

3.  05 

2.  97 

1.  12 

1.  21 

32 

3.  04 

2.  96 

1.  11 

1.  16 

44 

3.  04 

3.  00 

1.  10 

1.  12 

56 

3.  04 

3.  02 

1.  10 

1.  13 

68 

3.  04 

3.  04 

1.  09 

1.  11 

80 

3.  04 

3.  06 

1. 08 

1.  15 

92 

3.  04 

3.  01 

1.  07 

1.  14 

1504 

3.  04 

3.  10 

1.  07 

1.  15 

16 

3.  04 

2.  99 

1.  06 

1.  29 

28 

3.  04 

2.  86 

1.  05 

1. 23 

40 

3.  04 

2.  78 

1.  04 

1.  11 

52 

3.  04 

2.  79 

1.  04 

0.  97 

64 

3.  04 

2.  92 

1.  03 

0.  84 

76 

3.  03 

3.  04 

1.  03 

0.  88 

88 

3.  03 

3.  09 

1.  02 

0.  93 
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Table  II-l . (Continued) 


9 

i 


I 


/ -Iv 
(cm  ) 

'^smoothed 

^experimental 

^smoothed 

^experimental 

1600 

3.03 

3.  10 

1.01 

0.  96 

12 

3.  03 

3.  13 

1.  01 

0.  99 

24 

3.  03 

3.  15 

1. 00 

1.  04 

36 

3.  03 

3.  08 

1.  00 

1. 20 

48 

3.  03 

2.  90 

0.  99 

1.  04 

60 

3.  03 

3.  10 

0.  98 

0.  96 

72 

3.  03 

3.08 

0.  98 

1.  11 

84 

3.03 

2.  98 

0.  97 

1.  11 

96 

3.  03 

2.  92 

0.  97 

1.  05 

1708 

3.  03 

2.  94 

0.  96 

0.  98 

20 

3.  02 

2.  98 

0.  95 

0.  95 

32 

3.  02 

3.  03 

0.  95 

0.  96 

44 

3.  03 

3.  03 

0,  95 

1.  04 

56 

3.  02 

2.  95 

0.  94 

1.  01 

68 

3.  02 

2.  99 

0.  93 

0.  99 

80 

3.  02 

2.  88 

0.  93 

1. 04 

92 

3.  02 

2.  87 

0.  93 

0.  88 

1804 

3.  02 

2.  96 

0.  92 

0.  87 

16 

3.  02 

2.  98 

0.  91 

0.  89 

28 

3.  02 

2.  98 

0.  91 

0.  90 

40 

3.  02 

2.  98 

0.  90 

0.  90 

52 

3.  02 

2.  97 

0.  90 

0.  89 

64 

3.  02 

2.  97 

0.  89 

0.  88 

76 

3.  01 

3.  00 

0.  89 

0.  88 

88 

3.  01 

2.  96 

0.  89 

0.  90 

-40. 


I 


K 


1 


Table  II- 1.  (Concluded) 


V 

, -1,  snaoothed 

(cm  ) 

^^experimental 

smoothed 

^expe  rlment?  ’ 

1900 

3.01 

2.  96 

0.88 

0.84 

12 

3.  01 

2.  99 

0.  8 , 

0.84 

24 

3.01 

2.  98 

0.87 

0.85 

36 

3.01 

2.  95 

0.86 

0.83 

48 

3.  01 

3.  01 

0.86 

0.74 

60 

3.01 

3.  10 

0.  86 

0.84 

72 

3.01 

3.04 

0.85 

0.  85 

84 

3.01 

3.  07 

0.  85 

0.  85 

96 

3.01 

3. 05 

0.84 

0.  91 

2008 

3.01 

3.02 

0.84 

0.  91 

20 

3.  01 

3.  00 

0.  84 

0.  92 

32 

3.01 

2.  94 

6.  83 

0.  90 

44 

3.00 

2.  96 

0.83 

0.  85 

56 

3.00 

2.  97 

0.  82 

0,87 

68 

3.  00 

2.  94 

0.  82 

0.  90 

80 

3.  00 

2.86 

0.  82 

0.  85 

92 

3.  00 

2.  85 

0.  81 

0.75 

2104 

3.00 

2.  97 

0.81 

0.  66 

16 

3.  00 

3.03 

0.  80 

0.79 

28 

3.00 

2.  96 

0.80 

0.76 

40 

3.00 

3.  04 

0.80 

0.  76 

52 

3.00 

3.  00 

0.  79 

0.83 

64 

3.00 

2.  94 

0.79 

0.  81 

76 

3.00 

2.  92 

0.  79 

0.70 

88 

2.  99 

3.04 

0.  79 

0.74 

2200 

2.  99 

3.00 

0.  78 

0.  77 
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Table  II-2.  Reflectivity  Measurement  Errors  (rms) 


Wave  Number  Interval  AR,  Reflectivity  Error 


400  - 

600cm  ^ 

0.028 

600  - 

800cm"  ^ 

0.020 

800  - 

1000cm  ^ 

0.007 

1000  - 

1200cm  ^ 

0.  Oil 

1200  - 

1400cm  ^ 

0.011 

1400  - 

1600cm  ^ 

0.015 

1600  - 

1800cm"^ 

0.  013 

1800  - 

2000cm"^ 

0.  008 

2000  - 

2200cm"^ 

0.  01 1 

AR 


(n-8) 


There  is  no  a priori  reason  for  assurring  the  reflectivity  of 
carbon  should  be  smooth  as  a function  of  wave  numbei,  or  that  it  have 
the  particular  dependence  as  assumed  in  equation  II-7.  However,  the 
ARs  given  in  Table  II- 1 are  about  what  would  be  expected  when  detector 
sensitivity,  instrumental  response  and  atmospheric  absorption  are 
taken  into  consideration.  In  a sense  this  method  of  determining  reflec- 
tivity errors  is  conservative,  since  all  deviations  from  the  smoothed 
curve  are  treated  as  errors,  whereas  some  of  the  larger  deviations 
may  be  real. 


As  discussed  previously,  n and  k are  determined  from  the  gen- 
eralized Fresnel  equations. 


whe  re 


Q = 


(1  - Rj^)  cos  fi 
1 + Rj^-  2 y Rj^  cos 


P = 


- 2 /r  ^ cos  6 sin  4) 


1 + R 


,-z/r; 


cos  4: 


„2  ^2  2,2  . 2 „ , 

U - P = n - k - sin  6,  and 


QP  = nk, 

and  R is  the  reflectivity,  6 the  angle  of  incidence,  4^  the  phase  change 
on  reflection  (obtained  by  Kramers-Kronig  analysis),  and  n and  k the 
desired  optical  constants.  Q and  P,  and  therefore  n and  k,  are  deter- 
mined by  three  quantities;  Rj^,  4i  and  0,  and  therefore,  errors  in  these 
quantities  will  determine  the  probable  error  in  n and  k.  The  angle  of 
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Incidence,  6,  is  very  well  known  experimentally  and  analysis  has 
shown  that  the  error  in  the  measurement  of  fl(<  0.  2°)  is  negligible  in 
comparison  to  the  errors  in  R and  ip.  The  error  in  R,  AR,  was 
assumed  to  be 

AR  = ± 0. 015, 

and  is  representative  of  the  errors  given  in  Table  II-2.  The  errors  in 
ip  were  determined  by  examining  the  Kramers-Kronig  data  for  the 
experimental  and  smoothed  data,  in  a manner  analogous  to  the  analysis 
for  AR,  with  the  result: 

Aj/)  = 1.10  degrees. 

A computer  program  was  written  to  calculate  n and  k given  ip  and 
R,  and  then  four  calculations  were  carried  out,  setting 

R = R + AR,  and  ip  = ip  + Mp; 

R = R + AR,  and  ip  = ip  - Mp; 

R = R - AR,  and  ip  = ip  + l\ip; 

R = R - AR,  and  ip  = ip  ~ a<6; 

to  determine  n and  k.  The  two  extrema  (+  and  - about  the  input  value) 
were  then  selected  and  their  difference  from  the  input  value  was 
called  the  error  in  n and  k.  The  results  are  shown  in  Table  II-3,  for 
a range  of  n and  k from  zero  to  ten.  The  table  is  in  matrix  form, 
where  the  top  number  is  the  error  in  n and  the  bottom  number  is  the 
error  in  k.  As  an  example,  suppose  n was  measured  to  be  3.  0 and  k 
2.0,  then  Table  II-3,  gives  the  error  in  n as  ± 0.20  and  the  error  in 
k as  ± 0.  17.  Intermediate  cases  can  be  obtained  by  interpolation. 


! 


Table  II-3.  Errors  in  n and  k as  a Function  of  n and  k.  Assuming 


^ AR  = ±0.015,  Al/J  = ±1.1  , and  A9  = 0. 


0.0 

1 

1 

1.  0 

2.0 

3.  0 

5.0 

7.  0 

10.0 

An  = 

1.0 

± .05=:= 

. 09 

. 15 

.29 

.48 

. 85 

Ak  = 

± .05 

.09 

. 17 

.40 

.75 

1. 42 

Z.  0 

.07 

. 08 

. 12 

.21 

.42 

.66 

1.  12 

. 04 

.09 

. 10 

. 18 

.43 

. 80 

1.  58 

3.  0 

.11 

. 17 

. 20 

. 26 

. 52 

. 83 

1.  37 

.12 

. 17 

. 17 

. 17 

.44 

.83 

1.  64 

5.  0 

.22 

.34 

. 44 

. 53 

.73 

1.  10 

1.  81 

.36 

.40 

.44 

.44 

.40 

. 80 

1.  67 

7.  0 

.37 

. 55 

. 71 

. 86 

1.  12 

1.45 

2.  18 

.72 

. 74 

. 80 

. 83 

.78 

. 75 

1.  58 

10.  0 

.66 

. 92 

1.  17 

1. 41 

1.85 

2.26 

3.  06 

1. 52 

1.  55 

1.  53 

1.  59 

1. 62 

1.  54 

1.  53 

i j s *Top  figure  is  error  in  n,  lower  figure  is  error  in  k. 


Although  a definitive  analysis  was  not  performed,  the  calcula- 
tions from  which  Table  II-3  was  derived  indicated  that  the  errors  in 
n and  k were  more  sensitive  to  the  error  in  R as  opposed  to  Ip. 
Approximate  scaling  should  therefore  go  as  the  size  of  AR  (Table  II-2), 
although  it  should  be  kept  in  mind  that  the  Fresnel  equations  are  non- 
linear in  n and  k and  extreme  extrapolation  may  be  erroneous. 

The  optical  cr  stants  (n  and  k)  of  carbon  type  materials  (coals, 
ashes,  graphite,  etc.  ) are  the  subject  of  a vast  literature  and  the 
reader  is  referred  to  a review  by  Twitty  and  Weinman^^^^'  In 

order  to  illustrate  the  wide  variance  in  the  experimental  results,  the 
present  data  is  compared  with  three  previous  investigations  in  Figure 
II-9,  for  n t^'e  index  of  refraction,  and  in  Figure  II- 10,  for  k the 
extinction  coefficient.  In  each  figure  the  dashed  curves  are  the  results 
obtained  by  Foster  and  Howarth^^^^’  for  polycrystalline  graphite 

(method  of  preparation  not  discussed).  The  dotted  curves  are  the 
results  of  Lenham  and  Treherne^^^^'  for  pyrolytic  graphite 

O O 

(deposited  at  2280  C and  annealed  at  3200  C)  from  reflectivity  meas- 
urements with  the  electric  field  vector  polarized  perpendicular  to  the 
(presumably  oriented)  optic  axis  of  the  graphite  crystal(s).  The  solid 
curve  is  the  present  results  and  the  dot-dashed  curve  is  for  pressed 
carbon  soot  (obtained  for  an  oil  furnace  process,  fixed  carbon  per- 
centage = 98.5  percent,  particle  size  « 0.  023|om)  as  reported  by 
Foster  and  Howarth*-  ' •*.  The  results  shown  in  these  two 

figures  were  chosen  to  illustrate  the  wide  difference  in  results  obtained 
for  carbon  materials.  The  reader  is  referred  to  Ref.  11-18  for  a more 
complete  comparison.  Some  of  the  differences  are  probably  due  to 

experimental  error,  but  it  is  probable,  especially  from  the  investiga- 

r Ref  II-  1 9l 

tions  of  Foster  and  Howarth*^  ‘ that  most  of  the  apparent 

disagreement  is  real  and  due  to  such  factors  as  source  of  the  raw 
material,  and  method  of  preparation  of  the  samples. 
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2.  MgjO 

A single  crystal  of  MgO,  approximately  7.  6cm  long  by  3.  2cm 
wide,  was  wedged  to  prevent  internally  reflected  rays  from  entering 
the  interferometer,  and  polished  flat.  Reflectivity  measurements 

O O 

were  made  at  temperatures  of  300  K and  573  K.  The  reflectance 
spectra  are  shown  in  Figures  II- 11  and  11-12,  respectively.  The  peak 

o - 1 

in  the  300  K data  at  ~ 420cm  (one  data  point)  appears  to  be  spurious. 
The  reflectivity  curves  are  very  similar,  but  do  differ  in  two  details. 
First,  the  reflectance  of  the  hotter  case  is  in  general  lower  than  that 
for  the  cooler  case  in  the  region  of  400cm  ^ to  ~ 800cm  Second, 
the  slope  of  the  curve  at  about  v = 800cm  ^ is  somewhat  steeper  for 
the  cooler  measurement  compared  to  the  warmer  measurement. 

This  is  in  agreement  with  what  has  been  observed  previously. 

From  the  literature  and  our  data,  the  reflectance  changes  of  MgO  below 
400cm  ^ appear  to  be  significant,  and  indeed  the  infrared  allowed 
vibrational  mode  of  MgO  does  lie  below  400cm  Rather  than  extra- 
polate our  reflectance  data,  or  try  to  use  literature  reflectance  data 
for  the  wavenumber  region  below  400cm  ^ for  a Kramers-Kronig 
analysis,  we  analyzed  the  MgO  reflectance  data  using  classical  dis- 
persion analysis.  The  initial  values  of  the  three  dispersion  para- 
meters (Pj,  Vj  and  Yj)  were  taken  from  the  literature.  The  resulting 
best  fit  curves  to  the  input  reflectance  data  using  a single  oscillator 
are  shown  in  Figures  11-13  and  11-14,  and  are  to  be  compared  with 
Figures  II- 11  and  11-12,  respectively.  The  dispersion  parameters 
found  for  the  300°K  data  are  v = 389cm  \ Y = 71.3cm  ^ and 
4TTp  = 5.63.  The  parameters  for  the  573°K  data  are  v=  389cm  \ 

Y = 128.  2cm  ^ and  4np  = 5.  94.  For  both  temperatures  = 2.  3. 
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The  calculated  curves  are  very  smooth,  as  would  be  expected  for  a 
three  parameter  fit.  As  discussed  previously,  n and  k are  also 
determined  from  the  classical  dispersion  analysis,  and  they  are 
shown  in  Figures  11-15  and  II-16  for  T = 300  K,  and  in  Figures  11-17 
and  11-18  for  T = 573  °K.  In  Table  II-4,  the  index  of  refraction  (n) 
and  the  extinction  coefficient  (k)  of  MgO  are  tabulated  in  15cm 
increments  from  v = 400cm  ^ to  v = 2200cm  and  temperatures  of 
300°K  and  573°K. 


O 

The  present  results  (T  300  K)  are  compared  with  the  represen- 
tative data^^^^*  in  Figure  11-19.  Curves  A and  C are  the  present 

results  for  the  index  of  refraction  and  extinction  coefficient,  respec- 
tively. Curves  B and  D are  the  representative  data^^^^'  as 

O — X 

compiled  by  Whitson  (T  300  K).  Except  in  the  region  400cm  to 
500cm" ^ (a  region  of  low  s ignal-to-noise  for  our  apparatus).  The 
agreement  is  satisfactory  and,  in  fact,  is  typical  of  the  good  agree- 
ment between  previous  results  by  other  workers. 


3.  ZrOz 

The  reflectivity  spectra  of  ZrO^  between  400cm  ^ and  2200cm  \ 
at  a 45  degree  angle  of  incidence  were  measured  at  temperatures  of 
300°K  and  573°K.  These  reflectivity  spectra  are  shown  in  Figures 
11-20  and  11-21,  respectively.  A Kramers-Kronig  analysis  was  per- 
formed on  this  reflectivity  data  and  the  resulting  n and  k for  the  tern- 

o o 

peratures  of  300  K and  573  K are  shown  in  Figures  11-22  through  11-25. 

No  refractive  index  data  and  only  sparse  extinction  coefficient 
data  were  found  by  Whitson.  The  representative  reflec- 

O 

tance  data  given  by  Whitson  (T  = 300  K)  are  in  reasonable  agreement 
with  the  present  results.  The  representative  reflectance  data  of  Piriou 
[Ref.  11-22]  similar  spectral  structure  as  in  Figure  11-20,  but 

are  in  general  higher,  which  may  easily  be  attributed  to  a difference 
in  samples. 
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Figure  11-15.  The  Index  of 
dispersion  ai 


Ul: 


Table  II-4.  Index  of  Refraction  and  Extinction  Coefficient 

O O — 1 ~ 1 

for  MgO  at  300  K and  573  K from  400cm  to  2200cm 


T = 300  K 


T = 573  K 


V 

m"^) 

n 

k 

n 

k 

400 

3.  31 

4.  12 

2.88 

2.  97 

15 

2.35 

4.  09 

2.42 

3.04 

30 

1. 67 

3.  78 

2.  00 

2.  99 

45 

1.24 

3.  42 

1 . 66 

2.86 

60 

0.  96 

3.  09 

1.39 

2.  70 

75 

0.  78 

2.  79 

1.  18 

2.  52 

90 

0.  65 

2.  53 

1. 02 

2.  35 

505 

0.  55 

2.  31 

0.  90 

2.  18 

20 

0.48 

2.  11 

0.  80 

2.03 

35 

0.43 

1.  93 

0.  72 

1. 88 

50 

0.39 

1.  77 

0.  66 

1.74 

65 

0.  36 

1. 62 

0.  61 

1. 61 

80 

0.  33 

1. 48 

0.  57 

1.49 

95 

0.  31 

1. 35 

0.  54 

1. 37 

610 

0.30 

1.23 

0.  51 

1. 26 

25 

0.29 

1.  11 

0.49 

1.  15 

40 

0.29 

0.  99 

0.48 

1.05 

55 

0.29 

0.  88 

0.47 

0.  95 

70 

0.29 

0.  76 

0.  48 

0.  85 

85 

0.  31 

0.  64 

0.48 

0.  75 

700 

0.  34 

0.  53 

0.  50 

0.  66 

15 

0.  39 

0.42 

0.  52 

0.  57 

30 

0.46 

0.  32 

0.  56 

0.49 

45 

0.  54 

0.  25 

0.  60 

0.42 

60 

0.  61 

0.20 

0.  65 

0.  36 

1 


i 

> 

i ' 
r’' 


Table  II-4.  (Continued) 


T = 

300  °K 

T = 

573  °K 

V 

(cm  ) 

n 

k 

n 

k 

775 

0.  68 

0.  17 

0.  69 

0.  30 

90 

0.  74 

0.  14 

0.  74 

0.26 

805 

0.  79 

0.  12 

0.  78 

0.  23 

20 

0.  83 

0.  11 

0.  82 

0.  20 

35 

0.  88 

0.  10 

0.  86 

0.  18 

50 

0.  91 

0.  09 

0.  90 

0.  16 

65 

0.  95 

0.08 

0.  93 

0.  15 

80 

0.  98 

0.  07 

0.  96 

0.  13 

95 

1. 00 

0.  06 

0.99 

0.  12 

910 

1.  03 

0.  06 

1. 01 

0.  11 

25 

1. 05 

0.  05 

1. 03 

0.  10 

40 

1. 07 

0.05 

1. 05 

0.09 

55 

1. 09 

0.05 

1.  07 

0.  09 

70 

1.11 

0.  04 

1. 09 

0.  08 

85 

1.  13 

0.  04 

1.  11 

0.07 

1000 

1.  14 

0.04 

1.  13 

0.07 

15 

1.  16 

0.  03 

1.  14 

0.  07 

30 

1.  17 

0.  03 

1.  16 

0.  06 

45 

1.  18 

0.  03 

1.  17 

0.  06 

60 

1. 20 

0.  03 

1.  18 

0.05 

75 

1.21 

0.  03 

1.  19 

0.05 

90 

1. 22 

0.  03 

1. 20 

0.  05 

1105 

1.23 

0.  02 

1.21 

0.  05 

20 

1.24 

0.02 

1. 22 

0.  04 

35 

1. 25 

0.02 

1.23 

0.  04 

I ! 
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Table  II-4,  (Continued)  } 

i 

! 


T = 

300  °K 

T = 

573  °K 

V 

(cm“^) 

n 

k 

n 

k 

1150 

1. 26 

0.  02 

1.24 

0.  04 

65 

1.26 

0.  02 

1.25 

0.  04 

80 

1. 27 

0.02 

1.26 

0.  04 

95 

1. 28 

0.02 

1.27 

0.  03 

1210 

1. 29 

0.  02 

1. 27 

0.03 

1 

25 

1. 29 

0.02 

1. 28 

0.03 

i 

40 

1.  30 

0.  02 

1.29 

0.  03 

55 

1.  31 

0.  01 

1. 29 

0.  03 

70 

1.31 

0.  01 

1.  30 

0.03 

j 

85 

1. 32 

0.  01 

1.31 

0.02 

i 

1300 

1. 32 

0.01 

1.31 

0. 02  ; 

15 

1. 33 

0.  01 

1. 32 

0. 02  1 

30 

1. 33 

0.01 

1.  32 

0. 02  1 

J 

45 

1. 34 

0.  01 

1.33 

0. 02  j 

60 

1. 34 

0.  01 

1.  33 

0. 02  1 

i 

1 

1 

75 

1. 35 

0.  01 

1.34 

0. 02  1 

1 

90 

1. 35 

0.  01 

1.34 

0. 02  1 

1405 

1. 35 

0.  01 

1.  35 

0. 02  1 

! 

20 

1.  36 

0.  01 

1.  35 

0.  02 

35 

1. 36 

0.  01 

1.35 

0.  02 

50 

1.37 

0.  01 

1. 36 

0.  02 

■ 

65 

1. 37 

0.  01 

1.  36 

0.  02 

' 

80 

1. 37 

0.  01 

1.  36 

0.02  i 

i 

1 

95 

1. 38 

0.  01 

1. 37 

0.  02  ; 

1 

1 

1510 

1.  38 

0.01 

1. 37 

0.  02 

i 

Table  II-4.  (Continued) 


T = 

300  °K 

T = 

573  °K 

V 

cm  ^ ) 

n 

k 

n 

k 

1525 

1. 38 

0.01 

1. 37 

0.02 

40 

1.38 

0.01 

1. 38 

0.02 

55 

1.39 

0.01 

1. 38 

0.02 

70 

1. 39 

0.  01 

1. 38 

0.  02 

85 

1. 39 

0.  01 

1.39 

0.  01 

1600 

1.40 

0.  01 

1.39 

0.  01 

15 

1. 40 

0.01 

1. 39 

0.  01 

30 

1.40 

0.01 

1. 39 

0.01 

45 

1.40 

0.01 

1.40 

0.  01 

60 

1. 40 

0.01 

1. 40 

0.  01 

75 

1.41 

0.  01 

1.40 

0.  01 

90 

1.41 

0.  00 

1. 40 

0.  01 

1705 

1.41 

0.00 

1.41 

0.  01 

20 

1.41 

0.00 

1.41 

0.01 

35 

1. 42 

0.00 

1. 41 

0.  01 

1750 

1. 42 

0.00 

1.41 

0.  01 

65 

1. 42 

0.00 

1.41 

0.  01 

80 

1.42 

0.00 

1. 42 

0.  01 

95 

1. 42 

0.  00 

1.42 

0.01 

1810 

1. 42 

0.  00 

1. 42 

0.01 

25 

1.  43 

0.  00 

1.42 

0.  01 

40 

1. 43 

0.  00 

1. 42 

0.01 

55 

1.43 

0.00 

1.42 

0.  01 

70 

1.43 

0.00 

1.43 

0.  01 

85 

1.43 

0.  00 

1.43 

0.01 
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Table  II-4.  (Concluded) 


T = 

300°K 

T = 

573  °K 

V 

(cm  ) 

n 

k 

n 

k 

1900 

1.43 

0.  00 

1.43 

0.01 

15 

1.43 

0.  00 

1.43 

0.  01 

30 

1. 44 

0.  00 

1.43 

0.  01 

45 

1. 44 

0.  00 

1. 43 

0.  01 

60 

.1.44 

0.  00 

1.43 

0.  01 

75 

1. 44 

0.  00 

1. 44 

0.  01 

90 

1.44 

0.  00 

1.44 

0.  01 

2005 

1.44 

0.  00 

1. 44 

0.01 

20 

1.44 

0.  00 

1. 44 

0.01 

35 

1. 44 

0.  00 

1. 44 

0.  01 

50 

1. 45 

0.  00 

1. 44 

0.  00 

65 

1. 45 

0.  00 

1.44 

0.  00 

80 

1.45 

0.  00 

1. 44 

0.  00 

95 

1.45 

0.  00 

1.45 

0.  00 

2110 

1. 45 

0.  00 

1. 45 

0.  00 

25 

1. 45 

0.  00 

1.45 

0.  00 

40 

1.45 

0.  00 

1. 45 

0.  00 

55 

1.45 

0.  00 

1.45 

0.  00 

70 

1. 45 

0.  00 

1. 45 

0.  00 

85 

1. 45 

0.  00 

1. 45 

0.  00 

2200 

1. 46 

0.  00 

1.45 

0.  00 
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Figure  11-20.  The  reflectivity  spectrum  of  Zr02  li^  the  400cm  to  2000cm  region.  The  data  were 
obtained  at  an  angle  of  Incidence  of  45  degrees  and  at  a temperature  of  300  K.  Values 
belov/  400cm"  1 are  extrapolated. 


obtainsd  at  an  angle  of  incidence 
below  400cm- 1 are  extrapolated 


he  refractive  Index  (n)  of  Zr02  from  400cm  to  2000cm  as  determined  from  a 
ramer s-Kronlg  analysis  of  the  reflectivity  data  shown  in  Figure  11-21  (T  = 573  K 


A tabulation  of  n and  k for  Zr02  for  the  temperatures  T = 300°K 

O 

and  573  K is  given  in  Table  II-5  over  the  wavenumber  interval  from 
400cm  ^ to  2200cm  ^ in  15cm  ^ increments. 

4.  AI2O3 

The  reflectivity  spectrum  of  a single  crystal  sample  of  A1-0- 

- 1 -1  ^ ^ 

was  measured  in  the  wavenumber  region  400cm  to  2200cm  with 
the  electric  vector  polarized  perpendicular  to  the  C axis  of  the 
crystal  (EiC)  at  temperatures  of  300°K  and  Data  were 

obtained  also  for  E parallel  to  C,  but  appeared  to  be  overly  noisy  and 
were  not  analyzed.  The  reflectivity  spectra  for  T = 300 °K  and 
T = 678°K  are  shown  in  Figures  11-26  and  11-27,  respectively.  The 
spectra  are  very  similar,  but  in  general  the  higher  temperature 
sample  shows  lower  reflectivity  than  the  cooler  sample.  Also  for 
the  higher  temperature  sample  the  structure  in  the  reflectivity 
spectra  appear  broader  and  shifted  slightly  to  lower  wavenumber  in 
comparison  to  the  cooler  reflectivity  spectrum. 

A Kramers-Kronig  analysis  was  performed  on  each  reflectivity 
spectrum  and  the  derived  values  of  the  refractive  index  and  the 
extinction  coefficient  for  Al^O^  for  the  two  temperatures  are  shown 
in  Figures  11-28  through  11-31,  Table  II-6  lists  the  values  of  n and  k 
as  determined  in  this  investigation  of  AI2O2  for  the  two  temperatures 
in  the  wavenumber  interval  from  400cm"  ^ to  2200cm"  ^ in  15cm  '^ 
inc  rements. 

The  reflectivity  and  transmissivity  of  Al^O^  have  been  the  sub- 
ject of  many  investigations,  usually  with  the  object  of  determining  n 
and  k for  the  material.  Whitson^^^^’  has  examined  the  litera- 

ture and  made  a critical  selection  of  values  for  n and  k (for  nominally 

O 

300  K temperature).  These  are  shown  in  Figures  11-32  and  11-33,  and 
are  to  be  compared  with  Figures  11-28  and  11-29  (note  scale  difference), 
which  are  the  present  results.  These  values  of  n and  k were  then 
used  to  calculate  a reflectance  spectrum  for  Al^O^  at  an  incidence  angle 
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Table  II-5.  The  Refractive  Index  (n)  and  Extinction  Coefficient  (k) 
for  ZrO^  between  v = 400cm”  ^ and  2200cm  ^ for  T = 300  °K  and  573  °K. 


T = 300  K 


T = 573  K 


n 

k 

n 

k 

00 

2.  72 

2.  79 

1.28 

1. 76 

15 

1. 67 

2.77 

1. 26 

1. 57 

30 

1.48 

2.  07 

1.01 

1.26 

45 

1.  58 

1. 88 

1. 25 

1. 08 

60 

1. 62 

1.58 

1. 38 

1. 08 

75 

1.  94 

1. 59 

1. 42 

1.  16 

90 

1.  91 

2.  08 

1.29 

1.  18 

05 

1.39 

2.  01 

1. 22 

1.  17 

20 

1.24 

1. 84 

1.  14 

1. 06 

35 

1.  08 

1. 64 

1.  16 

0.  97 

50 

1.  11 

1.41 

1. 20 

0.  96 

65 

1.  38 

1. 37 

1.21 

0.  95 

80 

1. 22 

1. 74 

1.21 

1.  10 

95 

0.  83 

1.  55 

0.  95 

0.  90 

10 

0.  61 

1.34 

0.  97 

0.86 

25 

0.  69 

1.02 

0.  96 

0.  78 

40 

0.  78 

0.  95 

0.  96 

0.  75 

55 

0.  83 

0.  83 

0.  94 

0.  66 

70 

0.  92 

0.  78 

0.  92 

0.  61 

85 

0.  96 

0.  79 

0.  97 

0.  55 

00 

0.  96 

0.  77 

0.  98 

0.  54 

15 

0.  96 

0.  79 

0.  98 

0.  53 

30 

0.  86 

0.  79 

0.  93 

0.48 

45 

0.  79 

0.  69 

0.  94 

0.44 

60 

0.  77 

0.  58 

0.  91 

0.  38 
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Table  II- 5.  (Continued) 


i 

1 


I 


i 

! 


T = 

300  °K 

T = 

573  °K 

V 

(cm  S 

n 

k 

n 

k 

775 

0.  79 

0.48 

0.  95 

0.30 

90 

0.  82 

0.  39 

0.  99 

0.28 

805 

0.  84 

0.  31 

0.  99 

0.24 

20 

0.  93 

0.  22 

1. 04 

0.  19 

35 

1.03 

0.  19 

1.06 

0.  18 

50 

1.10 

0.  31 

1.  10 

0.  14 

65 

■ lvp3 

0.  18 

1.  15 

0.  16 

80 

1.  08' 

0.  16 

1.  15 

0.  16 

95 

1.  13 

0.  13 

1.  16 

0.  14 

910 

1.20 

0.  15 

1.  18 

0.  13 

25 

1.  19 

0.  15 

1.20 

0.  11 

40 

1. 22 

0.  12 

1.23 

0.  12 

55 

1. 26 

0.  13 

1. 22 

0.  13 

70 

1. 27 

0.  12 

1.21 

0.  10 

85 

1. 30 

0.  13 

1.26 

0.  06 

1000 

1.  32 

0.  13 

1. 30 

0.  08 

15 

1. 32 

0.  14 

1.31 

0.  11 

30 

1.  33 

0.  13 

1.  30 

0.  10 

45 

1. 34 

0.  14 

1. 33 

0.  09 

60 

1.  36 

0.  13 

1.32 

0.  11 

75 

1.  36 

0.  14 

1.34 

0.  11 

90 

1. 38 

0.  12 

1. 33 

0.  13 

1105 

1.39 

0.  13 

1.  34 

0.  10 

20 

1.41 

0.  14 

1. 37 

0.  12 

35 

1.42 

0.  14 

1. 36 

0.  15 

1 
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T = 

300°K 

T = 

573  °K 

r 

V 

(cm"^) 

n 

k 

n 

k 

• 

k 

1525 

1. 48 

0.  22 

1.38 

0.  15 

■ 

40 

1. 47 

0.  17 

1.33 

0.  07 

j:' 

i 

55 

1.  56 

0.21 

1. 48 

0.  03 

Y < 

70 

1.43 

0.  25 

1.40 

0.  16 

t 

V ' 

85 

1. 43 

0.  21 

1. 37 

0.  14 

[ ' 

1600 

1, 42 

0.  22 

1. 36 

0.  12 

k. 

15 

1. 40 

0.  19 

1. 38 

0.  11 

> i 

30 

1. 42 

0.  14 

1. 36 

0.  08 

> 

" » 

45 

1. 49 

0.  13 

1.  39 

0.  08 

f 

1 

60 

1. 45 

0.  21 

1. 45 

0.  08 

! 1 

75 

1. 42 

0.  15 

1. 43 

0.  19 

i 

90 

1.44 

0.  19 

1. 33 

0.  14 

) 

1705 

1.45 

0.  15 

1.41 

0.  05 

■! 

i 

1 

20 

1.47 

0.  20 

1. 42 

0.  16 

-i 

35 

1.42 

0.  16 

1. 35 

0.  10 

f 

50 

1. 44 

0.  15 

1. 39 

0.  09 

'] 

T- 

f 

65 

1.45 

0.  16 

1.40 

0.  12 

4. 

80 

1. 43 

0.  14 

1.37 

0.  12 

i 

95 

1. 45 

0.  14 

1.  38 

0.  10 

/ 

4 

J 

1810 

1.47 

0.  14 

1. 38 

0.  12 

t ' 

25 

1. 48 

0.  15 

1. 34 

0.08 

40 

1.48 

0.  18 

1. 38 

0.  06 

55 

1. 44 

0.  17 

1.41 

0.08 

1 

’:•  ; 

70 

1.47 

0.  14 

1.41 

0.  10 

! 

85 

1. 49 

0.  17 

1. 40 

0.09 
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Table  II-5.  (Concluded) 


V 

(cm"^) 


T = 300 °K 


T 


n 


k 


n 


1.46 

1. 45 

1. 46 
1.43 

1. 42 

1. 45 

1. 43 

1.47 

1. 48 
1. 45 

1. 42 

1. 45 

1. 43 

1. 46 

1. 46 

1. 45 

1. 44 

1. 44 
1. 42 

1. 45 


I /I  A 


t 


= 573 °K 


Figure  11-30.  The  refractive  index  (n)  of  AI2O3  in  the  region  400cm  to  2200cm"  as  determined  from  a 
Kramers-Kronig  analysis  of  the  reflectivity  data  shown  in  Figure  11-27  (T  = 678  K). 


Table  II-6. 

The  refractive 

index  (n)  and  extinction  co 

efficient  (k) 

r Al,Oo  in  the  region  400cm 

^ to  2200cm 

^ at  T = 300 

K and  678 

T = 300°K 

T = 678 

o 

K 

V 

(cm  ) 

n 

k 

n 

k 

400 

4.  70 

0.  00 

5.  86 

1. 32 

15 

4.  64 

0.  00 

4.  26 

1. 88 

30 

9.47 

0.  00 

6.  62 

2.  83 

45 

11.  59 

7.  77 

2.  46 

4.  50 

60 

1.  25 

6.  44 

1.  20 

2.41 

75 

1.  15 

2.  76 

1. 72 

1.  18 

90 

1.34 

1. 48 

2.  63 

0.  78 

505 

2.  66 

0.  70 

3.  02 

1. 01 

20 

3.40 

0.  71 

3.  27 

0.  87 

R • 

f , 

[ ' 

Table  II- 6.  (Continued) 


T = 

300  °K 

T = 

678  °K 

V 

(cm  ) 

n 

k 

n 

k 

775 

0.  08 

1.  67 

0.  30 

1. 45 

90 

0.  11 

1.  50 

0.  27 

1. 30 

805 

0.  12 

1. 39 

0.  28 

1.  17 

20 

0.  11 

1. 28 

0.  28 

1. 07 

35 

0.  10 

1.  13 

0.  26 

0.  94 

50 

0.  12 

1. 01 

0.  27 

0.  82 

65 

0.  11 

0.  89 

0.  27 

0.  68 

80 

0.  11 

0.  74 

0.  31 

0.  55 

95 

0.  13 

0.  58 

0.  36 

0.  42 

910 

0.  17 

0.  40 

0.45 

0.  29 

25 

0.  31 

0.  21 

0.  55 

0.  21 

40 

0.48 

0.  13 

0.  66 

0.  16 

55 

0.  62 

0.  10 

0.  75 

0.  14 

70 

0.  73 

0.  09 

0.  81 

0.  13 

85 

0.  81 

0.  09 

0.  86 

0.  10 

1000 

0.  88 

0.  10 

0.  92 

0.  10 

15 

0.  93 

0.  10 

0.  97 

0.  09 

30 

0.  96 

0.  11 

1. 00 

0.  09 

45 

1. 00 

0.  10 

1. 03 

0.09 

60 

1.  04 

0.  10 

1. 05 

0.  08 

75 

1. 07 

0.  10 

1.  08 

0.  07 

90 

1.  10 

0.  09 

1.  10 

0.  06 

1105 

1.  12 

0.  08 

1.  13 

0.  05 

20 

1.  16 

0.  08 

1.  16 

0.  05 

35 

1.  19 

0.  09 

1.19 

0.  05 
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Table  II-6.  (Continued) 


T = 

300  °K 

T = 

678  °K 

V 

(cm  ) 

n 

k 

n 

k 

1525 

1.  51 

0.  05 

1.49 

0.  00 

40 

1.  53 

0.02 

1.  52 

0.  00 

55 

1.  51 

0.  05 

1. 52 

0.  00 

70 

1.49 

0.  05 

1.  51 

0.  00 

85 

1.49 

0.  03 

1.  51 

0.  00 

1600 

1.49 

0.  01 

1. 51 

0.  00 

15 

1.  51 

0.  00 

1.  51 

0.  00 

30 

1.  54 

0.  00 

1.  53 

0.  00 

45 

1.  55 

0.  02 

1.  54 

0.  00 

60 

1.  54 

0.  02 

1.  54 

0.  00 

75 

1.  55 

0.  01 

1. 55 

0.  00 

90 

1.  56 

0.  01 

1.  55 

0.  00 

1705 

1.  57 

0.  03 

1.  56 

0.  00 

20 

1.  56 

0.  03 

1.  56 

0.  00 

35 

1.  55 

0.  04 

1.  55 

0.  00 

50 

1.  54 

0.  03 

1. 54 

0.  00 

65 

1.  54 

0.  00 

1.  53 

0.  00 

80 

1.  56 

0.00 

1.  55 

0.  00 

95 

1.  59 

0.  00 

1. 56 

0.  00 

1810 

1.  60 

0.  00 

1. 57 

0.  00 

25 

1.  60 

0.  02 

1.  58 

0.  00 

40 

1.  58 

0.  02 

1.  58 

0.  00 

55 

1.  58 

0.  00 

1.  58 

0.  00 

70 

1.  59 

0.  00 

1.  59 

0.  00 

85 

1.  60 

0.  00 

1. 61 

0.  00 
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Table  II-6.  (Concluded) 


T 

= 300 °K 

T = 

678  °K 

V 

(cm  ) 

n 

k 

n 

k 

1900 

1.  61 

0.  00 

1.  61 

0.  00 

15 

1. 62 

0.  01 

1.  61 

0.  00 

30 

1.  61 

0.  02 

1.  60 

0.  00 

45 

1.  59 

0.  00 

1.  60 

0.  00 

60 

1.  60 

0.  00 

1.  61 

0.  00 

75 

1.  61 

0.  00 

1.  62 

0.  00 

90 

1.  62 

0.  00 

1.  62 

0.  00 

2005 

1.  64 

0.  00 

1.  63 

0.  00 

20 

1.  63 

0.00 

1. 63 

0.  00 

35 

1.  63 

0.00 

1.  63 

0.  00 

50 

1. 64 

0.  00 

1.  62 

0.  00 

65 

1. 64 

0.  00 

1.  65 

0.  00 

80 

1.  65 

0.  00 

1. 66 

0.  00 

95 

1.66 

0.  00 

1.66 

0.  00 

2110 

1.  67 

0.  00 

1 . 66 

0.  00 

2125 

1.  67 

0.  00 

1. 66 

0.  00 

40 

1.  68 

0.  00 

1. 67 

0.  00 

55 

1.  69 

0.00 

1.  68 

0.00 

70 

1.  69 

0.  00 

1. 69 

0.  00 

85 

1.  70 

0.  00 

1.  69 

0.  00 

2200 

1.71 

0.  01 

1.71 

0.  00 

i 

i 

i 

s 

1 


J 
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by  Whitson  (Ref.  11-17)  from  the  literature. 


of  45  degrees.  The  result  is  shown  in  Figure  11-34,  which  is  to  be 
compared  with  Figure  11-26  (the  present  results).  The  agreement  is 
good  and  perhaps  fortuitous  considering  the  fact  that  samples  obtained 
from  different  sources,  or  prepared  differently  often  show  significant 
differences  in  their  reflectivity  spectra. 


In  an  extensive  investigation  of  the  room  temperature  (=“  300  K) 
reflectivity  spectrum  Barker^ examined  eight  different  sam- 
ples obtained  from  various  sources.  The  reflectivity  data  were 
reduced  using  a classical  dispersion  analysis.  Chromium  impurities 
even  up  to  concentrations  of  2.4  X 10^^  ions  per  cc  were  not  detected 
in  the  infrared  spectrum,  however,  the  method  of  preparation  of  the 
reflecting  surfaces  did  have  a significant  effect  on  the  reflectivity 
spectrum,  and  is  illustrated  by  a quote  from  this  investigation^ 


"An  etch  of  molten  boron  oxide  and  lead  oxide  was  used  to 
remove  a few  microns  from  the  surfaces  of  samples  3 and  5. 
Spectacular  changes  (in  the  reflectivity  spectrum)^!-  resulted. 
Etching  the  surface  removes  nearly  all  traces  of  the  forbidden 
modes.  Grinding  with  diamond  dust  by  the  usual  techniques  then 
caused  the  modes  to  reappear.!  Grinding  with  course  grit,  15|j,ni 
particle  size,  caused  the  forbidden  modes  to  reappear  with 
larger  strength.  Further  grinding  with  6|j,m  grit  weakened  the 
forbidden  modes.  " 


It  would  thus  appear  (at  least  for  Al^O^)  the  method  of  surface 
preparation  has  a larger  effect  than  impurity  concentration  on  the 
reflectivity  spectrum.  Our  sample  was  polished  with  a very  fine  grit 
(~0.3um)  and  in  the  course  of  the  experiments  heated  up  to  about 
700 °K  which  should  have  relieved  to  a large  degree  any  of  the  strain 
which  allow  the  appearance  of  the  forbidden  modes. 


=i=Words  in  parenthesis  added  for  clarification  in  quotation. 


Figure  11-34.  Calculated  reflectivity  at  45  degrees  incidence  for  AI2O3  (T  300  K)  in  the  region 

100cm"  1 to  2000cm"  1.  The  n and  k values  used  are  the  representative  values  from  the 
compilation  by  Whitson  (Ref.  11-17). 


D. 


SUMMARY  AND  GENERAL  CONCLUSIONS 

The  reflectLvlty  spectra  of  Al^O^,  MgO  and  ZrO^  at  room  tem- 
perature and  elevated  temperatures  (~  600  K),  and  carbon  at  room 
temperature  have  been  measured.  From  this  data  the  real  and 
Imaginary  parts  of  the  Index  of  refraction  have  been  determined  by  a 
Kramers-Kronlg  analysis,  or  classical  dispersion  theory.  Agreement 
with  previous  (mostly  room  temperature  data)  Investigations  Is  good. 

An  error  analysis  based  on  the  carbon  data  has  been  carried  out.  It 
Is  concluded  that  the  reflectance  method  for  determining  n and  k Is 
fairly  accurate  as  long  as  n and/or  k are  not  too  large  or  too  small. 

In  particular,  the  method  Is  very  poor  lor  determining  k ^ 0.05. 

Direct  determination  of  n and  k from  the  Fresnel  equation  was 
found  to  be  Inadequate  for  many  situations.  Small  changes  In  the 
reflectivity  data  often  produced  large  changes  In  the  deduced  values 
of  n and  k.  Kramers-Kronlg  analysis  provided  a much  more  stable 
solution  to  the  problem. 

Reflectivity  data  at  higher  temperatures  were  not  obtained,  due 
to  heater  and  sample  failure.  This  problem  Is  not  considered  serious 

O 

and  the  higher  temperature  data  1000  K)  should  be  obtainable. 

As  discussed  above,  the  straightforward  reflectance  type  meas- 
urements described  In  this  Investigation  give  very  poor  accuracy  for 
k < 0.  05.  However,  by  coating  the  sample  to  be  measured  with  a 
dielectric  of  suitable  Index  of  refraction,  Fahrenfort^^®^'  has 

described  and  developed  a technique  called  "Attenuated  Total  Reflec- 

_4 

tlon"  which  appears  to  be  capable  of  measuring  k down  to  ~ 10 
Our  apparatus  appears  to  be  suitable  for  making  such  measurements, 
but  would  require  some  exploratory  and  developmental  effort  to  verify 
the  validity  of  the  technique  (especially  for  elevated  temperature  meas- 
urements). Small  values  of  k can  also  be  determined  by  measuring 
transmissivity  spectra.  This  technique  has  been  well  developed  and 
xi.easurements  at  elevated  temperatures  have  been  successfully 
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SECTION  III 


THEORETICAL  COMPUTATION  OF  PARTICLE  EMISSIVITY 


One  requirement  of  the  Particle  Optical  Properties  Measurements 
(POPM)  Program  is  to  calculate  the  emissivity  in  the  infrared  region 
of  small  spherical  particles  from  the  complex  index  of  refraction  of 
the  material  composing  the  particle.  It  is  the  purpose  of  this  section 
to  outline  the  derivation  of  mathematical  expressions  necessary  to 
make  these  calculations  and  to  apply  them  to  predict  the  emissivity 
for  particles  of  the  materials  for  which  bulk  refractive  indices  are 
presented  in  Section  II. 


A.  DERIVATION  OF  THEORETICAL  FORMULAS  FOR  PARTICLE 
EMISSIVITY 

At  temperatures  above  absolute  zero,  all  matter  radiates  and 
absorbs  energy  in  the  form  of  electromagnetic  waves  distributed  over 
a spectrum  of  frequencies.  The  well  known^^®^’  theoretical 

expression  for  the  power  radiated  from  the  surface  of  a radiating 
body  is; 

R(v,  T)  dv  = e(v,T)  (2rrhv^/c^)  [exp(hv/kT)  - l]  ^ dv  (III-l) 


where 

k 

h 

V 

c 

T 

€(v,  T) 


= Boltzmann's  constant, 
= Planck's  constant, 

= frequency, 

= velocity  of  light, 

= temperature,  and 
= emissivity. 
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and  R(v,  T)  dv  Is  the  power  per  unit  area  emitted  Ln  the  frequency 
range  from  v to  v + dv.  Multiplying  this  by  the  surface  area  of  a 
sphere  of  radius  r,  we  obtain  the  total  radiated  power  emitted  by  the 
sphere  in  the  frequency  range  v to  v + dv, 

P(v,  T)dv  = e(v,  T)  Trr^{8rrhv^/c^)  [exp(hv/kX)  - 1]"^  dv.  (III-2) 

The  emissivity  appears  as  a factor  in  this  expression.  Thus,  equation 
III-Z  can  be  regarded  as  the  basic  definition  of  the  emissivity  of  a 
spherical  particle  of  radius  r. 

The  spectral  distribution  of  radiated  power,  P(v,T),  can  be 
computed  theoretically  by  applying  Kirchhoff's  law  and  solving 
Maxwell's  equations.  Kirchhoff's  law  states  that  under  conditions  of 
thermal  equilibrium,  the  power  emitted  in  any  small  frequency  range 
is  equal  to  the  power  absorbed  in  the  same  frequency  range.  The 
power  absorbed  is  then  computed  by  solving  Maxwell's  equations  for  the 
scattering  of  electromagnetic  waves  from  a spherical  particle.  By 
comparing  the  expression  we  obtain  by  this  method  with  equation  III-2 
we  obtain  a theoretical  formula  for  the  emissivity.  A recent  detailed 
study  by  Baltes^^^^’  has  shown  the  emissivity  obtained  by  these 

thermal  equilibrium  arguments  is  also  appropriate  for  conditions 
when  the  radiating  object  and  the  radiation  field  are  not  in  thermal 
equilibrium.  In  the  remaining  part  of  this  section  we  will  carry  out 
the  procedure  outlined  above  to  obtain  formulae  for  computing  emis- 
sivity. 

Imagine  a spherical  particle  suspended  inside  a hollow  enclosure 
(hohlraum).  The  entire  system  which  consists  of  the  walls  of  the 
enclosure,  the  particle  and  the  radiation  field  inside  the  enclosure  is 
assumed  to  be  in  thermal  equilibrium  at  a temperature  T.  The  volume 
density  of  radiant  energy  between  the  frequencies  v and  v + dv  is  given 
by  the  Planck  formula, 


r 


i 
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(III-3) 


U(v,  T)  dv  = (8nhv^/c^)  [exp(hv/kT)  - 1]'^  dv. 

This  radiation  field  is  homogeneous  and  isotropic.  It  can  be  analyzed 
into  a sum  of  individual  unpolarized  plane  waves  traveling  with  the 
speed  of  light,  c,  in  all  possible  directions.  Let  us  establish  a coord- 
inate system  with  its  origin  at  the  center  of  the  particle.  The  intensity 
of  the  radiation  traveling  in  the  direction  of  a unit  vector  specified  by 
the  polar  coordinates  6 and  cp  and  contained  within  the  infinitesimal 
solid  angle  df2  is, 

dl(v,  T)  = K(v,T)dO.  (III-4) 

The  quantity  K(v,  T),  defined  by  equation  III-4,  has  no  angular  depend- 
ence and  can  be  related  to  the  energy  density  by  the  following  analysis. 

The  intensity  of  a plane  wave  is  related  to  the  energy  density  in  that 
wave  by  the  expression. 


I(v,  T)  = cU(v,  T).  (III-5) 

Combining  equations  III-4  and  III-5  and  integrating  over  all  solid  angles 
gives  the  desired  relationship  between  K(v,  T)  and  U{v,  T), 

K(v,  T)  = (c/4tt)  U{v,  T).  (III-6) 

Each  plane  wave  of  intensity  dl(v,  T)  is  scattered  and  partially 
absorbed  by  the  spherical  particle.  This  process  is  most  conveniently 
expressed  in  terms  of  a cross  section.  The  amount  of  power  in  the 
range  v to  v + dv  absorbed  by  the  particle  from  a plane  wave  of  inten- 
sity dl(v,  T)  is  given  by, 

dPabs(\>.  T)  dv  = a(v,  T)  dl(v,  T)  dv,  (III-7) 

where  ct(v,  T)  is  the  absorption  cross  section.  Substitute  equation  III-4 
into  equation  III-7  and  integrating  over  all  solid  angles, 


(III-8) 


Pabs^^’ = 4na(v,  T)  K(v,  T)  dv. 

Then  substitute  equation  III-6  Into  equation  III-8, 

Pabs^^’ = ca(v,  T)  U(v,  T)  dv.  (HI-9) 

Utilizing  Klrchhoff's  law  which  equates  the  power  absorbed, 

Pabs^'^’  power  emitted,  P(v,  T)  dv,  and  substituting  from 

equation  III-3  we  obtain, 

P(v,  T)  dv  = a(v,  T)  (8nhv^/c^)  [exp(hv/kT)  - l]~^  dv.  (III-IO) 

Comparing  equations  III- 10  and  III-2  we  see  that  the  emissivity  is 
given  by. 


I 


e{v,T)  = a(v,  T)/rrr^,  (III-ll) 

i.e.,  it  is  just  the  absorption  cross  section  divided  by  the  geometric 
cross  section.  Van  de  HulsJ^®^'  calls  this  ratio  of  cross 

sections  the  absorption  efficiency  and  has  denoted  it  by  the  symbol 
Qabs‘  Thus,  the  emissivity  of  a spherical  particle  is  equal  to  its 
absorption  efficiency, 

e(v,T)  = (III-12) 

Van  de  Hulst^^^^'  gives  the  following  set  of  Mie  theory 

formulas  for  the  calculation  of  Q , . 

abs 

00 

^abs  = ^ (2n  + 1)  [Re(a^  + bj  - |a^|2  . \hj],  (HI-13) 

n= 1 
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where 


X 

= Znr/X, 

(Ill- 14) 

X 

= c/v  = wavelength. 

(III- 15) 

a 

n 

jn(y)  in(x)  - mj^(y)  j^(x) 
jn(y)  (x)  - mj^(y)  h^^^  (x) 

(III- 16) 

b 

n 

rnjn(y)  jjx)  - j^(y)  j|^(x) 
injn(y)  hjf  ^x)  - j^(y)  hj^^^  (x) 

(III- 17) 

m 

= n(v,  T)  - ik(v,  T), 

(III- 18) 

n(v,T) 

= index  of  refraction. 

k(v,  T) 

= coefficient  of  extinction,  and 

y 

= mx. 

(III- 19) 

and  h|^^^(x)  are  the  Ricatti-Bessel  functions.  The 

primes  in  equations  III- 16  and  III- 17  denote  differentiation  with 
respect  to  the  argument  of  the  function.  The  temperature  and  fre- 
quency dependence  of  [and  thus  of  e(v,  T)]  is  due  to  the  depend- 

ence of  the  index  of  refraction  and  the  coefficient  of  extinction  on 

these  quantities. 


B.  . CALCULATIONS  OF  THE  EMISSIVITY  OF  SMALL  SPHERICAL 
PARTICLES 

A computer  program  entitled  "SPECTRUM”,  which  calculates 
the  emissivity  by  the  method  above  has  been  written.  The  input  to 
the  program  is  the  refractive  index,  n(v,  T),  the  extinction  coefficient, 
k(v,  T),  and  the  radius  of  the  particle.  The  program  computes  and 
plots  the  emissivity  as  a function  of  wavenumber  (cm'^). 


-101- 

— 

iaiaMiii'ifiiiffii  I 111  t iniTiiiiii  ■ 


I 


r 


k 


I 


Shown  in  Figures  III- 1 through  IIl-lO  are  the  calculated  emissivi- 
ties  of  spherical  carbon  particles  whose  radii  vary  from  0.  5|im  to  50pm. 

The  wavenumber  range  of  the  calculations  is  300  to  1800cm  . The  n 

and  k values  used  are  from  the  smoothed  data  for  carbon  (Section  II-C.  1) 
except  for  Figure  III-8,  which  used  the  n and  k derived  from  the  experi- 
mental data.  This  figure  is  to  be  compared  with  Figure  III-7.  The  dif- 
ferences between  the  two  curves  are  noticeable,  but  quite  small. 

As  was  expected,  the  variation  of  emissivity  with  wavenumber 
is  very  smooth  and  rather  slowly  varying.  The  change  in  emissivity 
with  particle  size  is  much  more  dramatic,  and  in  fact,  for 
1.0p,m  s r ^ lOpim  sized  particles  the  calculated  emissivity  is  greater 
than  one  over  a large  portion  of  the  wavenumber  interval  plotted.  At 
first  glance  this  appears  surprising,  however,  it  should  be  remem- 
bered that  these  sized  particles  are  in  the  "Mie"  region,  i.  e. , their 
circumference  is  of  the  order  of  magnitude  of  the  wavelength  of 
radiation.  As  is  well  known,  particles  in  this  size  region  frequently 
have  scattering  cross  sections  larger  than  their  geometrical  cross 
sections,  due  to  resonance  between  vibration  modes  of  the  particle 
and  the  radiation  and  thus,  absorption  or  em/ission  cross  sections 
larger  than  the  geometric  cross  section  shc.uld  not  be  too  surprising. 

As  the  particle  size  becomes  larger  (i.  g.  , Figures  III-9  and  IIl-lO), 
the  calculated  emissivity  falls  below  one  and  is  more  or  less  constant 
above  500cm  ^ (shorter  than  X = 20p,m).  Thus,  the  larger  carbon 
particles  act  like  classical  greybodies,  as  they  should. 

For  particles  large  compared  with  the  wavelength  of  radiation, 
the  computed  particle  emissivity  should  approach  the  emissivity  of  a 
plane  bulk  sample.  The  emissivity  at  normal  incidence,  e^,  is  par- 
ticularly easy  to  compute  for  an  opaque  material  such  as  carbon. 

= 4n/[(n  + 1)^  + k^]  (111-20) 
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Figure  III-2.  Calculated  emisslvlty  of  a spherical  carbon  particle,  radius  = l.Opm. 
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Calculated  emissivity  of  a spherical  carbon  particle,  radius  = 20pm 
(Experimental  refractive  index  data.  ' 


Calculated  emlssivlty  of  a spherical  carbon  particle,  radius  = SOiJjn 
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With  n = 3 and  k = 1,  then  e = . 71  in  reasonable  agreement  with  the  i 

' n ® I 

limiting  values  of  carbon  particle  emissivity  computed  at  higher  fre-  ^ 

5 

quencies  for  the  larger  particles.  { 

i 

Shown  in  Figures  III- 11  through  III-28  are  the  calculated  emis-  I 

sivities  for  spherical  Al^O^  particles  of  radii  from  0.  5gm  to  50gm.  I 

Two  sets  of  results  are  shown.  The  upper  figure  on  each  page  gives  *| 

the  emissivity  as  calculated  from  literature  values  of  n and  k (i.e., 

the  representative  values  from  the  handbook  compiled  by  Whitson,  f 

Ref.  III-5).  The  lower  curves  were  calculated  from  the  n and  k 5 

determined  in  this  investigation.  The  calculated  emissivity  spectra  i| 

of  spherical  Al^O^  particles  are  very  structured  in  the  400  to  900cm 

region.  =!=  Part  of  this  structure  is  due  to  the  structure  of  the  functions  ^ 

n(v,  T)  and  k(v,  T),  however,  for  very  small  particles  some  of  the  I 

structure  is  also  a manifestation  of  the  polariion  modes^^®^'  of  i 


the  particle.  The  calculated  emissivities  are  seen  to  be  in  good  \ 

agreement  for  v less  than  about  1000cm  In  this  region  the  emis- 

sivity  spectra  calculated  from  the  literature  n and  k values  appear  i 

sharper,  and  in  some  regions  more  intense  (i.e.,  large  emissivity)  3 

? 

than  that  calculated  from  the  n and  k obtained  in  the  present  investi-  } 

gation.  This  may  well  be  due  to  a resolution  effect.  In  any  case  ’ ; 

the  behavior  of  the  emissivity  in  the  wavenumber  region  400cm  ^ it 

(25ym)  to  1000cm  ^ (lOgm)  is  seen  to  depend  critically  on  the  radius 

of  the  emitting  particle.  J 

; j 

*In  the  region  v = 400cm  ^ and  lower,  our  results  depend  on  how  our  i 

reflectivity  spectra  were  extrapolated  for  v s 400cm"l,  and  thus  ; I 

comparison  of  the  two  calculations  is  not  valid  in  this  region.  > 

_ I ' 

**The  resolution  employed  in  this  investigation  was  ~ 15cm  . The 
resolution  applicable  to  the  literature  values  of  n and  k is  not 

known.  If  the  resolution  for  the  literature  values  of  n and  k is  1 

better  (i.e.,  < 15cm"^)  then  the  more  complex  structure  seen  in 

the  upper  curves  is  probably  real,  however,  if  the  resolution  for 

the  literature  values  of  xi  and  k is  poorer,  then  the  more  complex 

structure  may  be  due  to  an  artifact,  such  as  the  way  the  data  for 

n and  k were  digitized  for  compilation  in  the  handbook  (Ref.  Ill- 5) 

by  Whitson. 
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Figure  III- 18.  Calculated  emissivLty^for  a 3(j,m  radius  Al^O-  particle,  n and  k present  results. 
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The  comparison  of  the  calculated  particle  emissivities  for 
v>  lOOOcm'^  is  not  very  satisfactory,  and  in  fact,  becomes  progres- 
sively poorer  as  the  particle  radius  increases.  This  is  due  almost 
solely  to  the  differences  in  the  k (extinction  coefficient)  values  as 


compiled  in  the  handbook  by  Whitson  and  the  present  results  (see 
Section  II,  Table  II-6).  At  about  1000cm"  ^ (10vj,m)  the  literature 
value  of  k is  ~ 3 X lO"^,  decreasing  monotonically  to  ~ 3 x 10 


2000cm"  (5vi,m).  The  present  experimental  results  give  k in  the 

range  1 X lO"^  to  5 X lO"^  for  lOOOcm"^  s v ^ 2000cm"  except  for 


a 


few  values. 


Theoretically,  for  a perfect  crystal,  for  frequencies  greater 
than  the  highest  infrared  allowed  lattice  resonance  frequency  the 
extinction  coefficient  should  become  small  until  electronic  effects 
become  important.  In  fact,  had  we  made  a classical  dispersion  fit  to 
our  Al^Oj  reflectivity  data  using  literature  or  near  literature  values 
for  the  dispersion  parameters,  we  would  have  obtained  considerably 
lower  values  of  k for  v>  1000cm" From  a consideration  of  the 
noise  encountered  in  our  experimental  reflectivity  spectra  (see 
Section  II-C.  1 on  carbon,  where  the  effect  of  experimental  error  on 
the  determination  of  n and  k is  discussed)  we  conclude  errors  of 
~ 0.  05  in  k are  quite  possible.  We  feel  that  the  single  crystal  of 

AFO-  used  in  this  investigation  was  sufficiently  pure,  and  lacking  in 
2 3 _ 2 

imperfection,  so  that,  in  the  region  v>  1000cm  , k should  be  con- 

siderably smaller  than  the  measured  value.  Thus,  the  calculated 
particle  emissivities  in  this  region  should  be  considered  upper  limits. 

For  an  imperfect  crystal,  however,  the  conclusions  can  be 
quite  different.  The  presence  of  impurity  sites,  vacancies,  disloca- 
tions, etc.  , can  in  general  break  down  the  selection  rules  applicable 
to  a perfect  crystal.  Thus,  resonances  not  previously  allowed  to  be 
infrared  active  could  become  infrared  active.  In  such  cases  the  value 
of  k could  be  significantly  perturbed  upward,  and  the  emissivity  of  an 
spherical  imperfect  crystalline  particle  would  then  exhibit  much  higher 
emissivity. 


1 
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Figure  III-29  is  a plot  of  the  emissivity  of  Al^O^  which  was 
calculated  using  the  simple  formula  shown  in  equation  111-20,  and  the 
literature  values  for  n and  k.  This  plot  should  be  compared  to 
Figure  III-27  which  was  calculated  using  exact  Mie  theory  for  a 
particle  of  radius  50y,m.  It  is  apparent  that  this  simple  limiting 
formula  has  reproduced  the  complicated  shape  of  the  emissivity 
quite  well  in  the  region  of  frequency  less  than  1000cm  The  only 
feature  it  failed  to  reproduce  is  the  sharp  minimum  at  about  350cm  , 

however,  this  is  at  a wavelength  of  the  same  order  as  the  particle 
size  and  so  the  simple  formula  is  not  expected  to  be  valid.  Mie  theory 
will  predict  zero  emissivity  if  the  extinction  coefficient,  k,  is  zero 
whereas  the  simple  formula  can  predict  a finite  emissivity.  This 
will  only  cause  difficulty  when  k is  zero  or  extremely  small  and  is 
the  reason  the  emissivity  falls  off  at  high  frequency  in  Figure  III-27 
and  does  not  in  Figure  III-29.  Keeping  in  mind  its  limitations, 
equation  111-20  is  an  extremely  useful  formula,  adding  physical  insight 
and  enabling  one  to  do  "back  of  the  envelope"  calculations. 

For  completeness,  particle  emissivity  calculations  for  spherical 
particles  (0.  5p,m  ^ r ^ SO^im  range)  were  carried  out  for  Al^O^ 

(T  = 678°K),  MgO  (T  = 300 °K  and  573°K),  and  2.rO^  (T  = 300"k  and 
573°K). 

The  results  for  Al^O^  (T  = 678  °K)  are  shown  in  Figures  111-30 
through  III-38.  The  results  are  seen  to  be  quite  different  than  obtained 
from  the  T = 300 °K  data,  and  reflect  the  changes  in  n and  k deduced 

O . O 

from  the  reflectivity  spectrum  measured  at  T = 300  K and  678  K. 

The  calculated  emissivities  for  spherical  particles  (0.  5p,m  s r 
^ 50vurv  range)  of  MgO  are  shown  in  Figures  III-39  through  111-47  for 
T = 300'k,  and  Figures  111-48  through  III-56  for  T = 573  K.  The  MgO 
emissivity  spectra  are  much  simpler  than  those  calculated  for  Al202» 
This  is  because  the  MgO  cfystal  structure  is  much  simpler  than  the 
AI2O2  crystal  structure. 
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Figure  III- 30.  Calculated  emlsslvlty  for  a 0.  5gjn  radius  A1~0_  particle,  T = r>78 


Calculated  emisslvity  for  a 50.  Oy,m  radius  Al;>Oo  particle 


Calculated  emisslvity  for  a l.Opjn  radius  MgO  particle,  T = 300  K 


Calculated  emissLvity  for  a 3.  Op,m  radius  MgO  particle,  T = 300  K 


Calculated  emissivlty  for  a 30.  Ojim  radius  MgO  particle,  T = 300  K 


Calculated  emlsslvlty  for  a 1.0pm  radius  MgO  particle,  T = 573  K 


Figure  III-50.  Calculated  emisslvity  for  a 2.  0|jjn  radius  MgO  particle,  T = 573 


venumbe 


Figure  III-55.  Calculated  emissivity  for  a 30.  Op.m  radius  MgO  particle,  T - 573 


The  calculated  emissLvitles  for  spherical  particles  (0.  5u,m  s r 
s 50[jjn)  for  ZrO^  are  shown  In  Figures  III-57  through  III-65  for 
T = 300°K,  and  Figures  III-66  through  III-74  for  T = 573  °K.  Since 
the  reflectivity  data  were  obtained  for  a polycrystalline  sample,  the 
results  are  applicable  to  spherical  polycrystalline  particles  of  ZrO^ 
and  not  the  single  crystalline  form. 


C.  SUMMARY  AND  GENERAL  CONCLUSIONS 

The  theory  of  small  particle  emissivity  has  been  summarized. 
The  emissivity  of  a spherical  particle  is  equal  to  its  absorption 
efficiency  as  calculated  from  Mie  theory.  A computer  program  has 
been  written  which  calculates  the  emissivity  of  spherical  particles. 
Using  the  data  given  in  Section  II,  spherical  particle  emissivities  have 
been  calculated  for  particle  sizes  in  the  range  0.  5(j,m  to  50p,m.  For 
small  particles  (<  10p,m)  the  emissivity  often  exceeds  one,  since  the 
absorption  cross  section  is  larger  than  the  geometrical  cross  section. 
For  larger  particles  the  calculated  emissivity  does  indeed  reduce  to 
its  classical  value.  The  emissivity  of  spherical  particles  is  found  to 
be  very  dependent  on  the  values  of  n and  k,  as  well  as  the  size  (dia- 
meter) of  the  particle.  Also,  the  calculated  emissivity  of  intermed- 
iate (lOgm  s r ^ 50p,m)  sized  particles  becomes  quite  large  (ca  1)  for 
relatively  small  values  of  k (~0.02),  which  in  general  are  poorly 
determined  by  reflectance  type  measurements.  Furthermore,  these 
values  of  k which  can  lead  to  high  emissivity  values  under  proper 
conditions  are  much  more  likely  to  occur  for  the  impure  and  imper- 
fect crystals  in  particular  than  for  large,  carefully  prepared  bulk 
samples  where  k < 0.001  at  frequencies  greater  than  the  last  infrared 
allowed  frequency. 

The  calculated  emissivity  of  small  spherical  particles  is  seen 
to  be  critically  dependent  on  size,  and  complex  index  of  refraction 
(i.  e.,  n and  k).  We  therefore  conclude  that  real  particles,  even  if 
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they  are  quasi- spherical,  will  probably  have  an  emissivity  quite  dif- 
ferent from  that  calculated  from  Mie  theory,  since  the  bulk  values  of 
n and  k may  not  be  applicable  to  very  small  particles.  Measurement 
of  individual  particle  emissivities  are  still  required. 


Figure  III-58.  Calculated  emissivity  for  a 1.0|jm  radius  ZrO-,  particle,  T = 300 


ulated  emLssivLty  for  a 2.  Oyrn  radius  ZrO^  particle,  T = 300  K 


Calculated  emissivity  for  a 3.  Ogrn  radius  ZrO  particle,  T - 300  K 


Wavenumber  (cm~  ) 

Figure  III-6Z.  Calculated  emissLvlty  for  a 10.  Oum  radius  ZrO-,  particle,  T = 300 


Figure  III-68.  Calculated  emissivLty  for  a 2.  Oum  radius  ZrO  particle,  T = 573 
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SECTION  IV 


! 


SINGLE  PARTICLE  EMISSIVITY  EXPERIMENT  if 


A.  EXPERIMENTAL  APPROACH 

The  basic  objective  of  the  POPM  program  is  the  comparison  of 
particle  emissivities  computed  from  bulk  data  with  values  measured 
for  single  particles.  This  section  describes  the  progress  which  has 
been  made  experimentally  in  holding  and  heating  a single  particle, 
then  observing  the  infrared  radiation  from  the  hot  particle,  and 
finally  analyzing  the  infrared  data  to  obtain  particle  emissivities. 

The  single  particle  emissivity  experiment  is  based  on  simple 
scientific  principles  which  require  a combination  of  state-of-the-art 
techniques.  A single  particle  of  Al^O^,  MgO,  ZrO^  or  carbon  is 
electrostatically  suspended  at  the  focal  point  of  a Nd:YAG  laser,  which 
is  used  for  heating  the  particle.  Conventional  optics  are  used  to 
collect  and  spectrally  analyze  the  radiation  from  the  heated  particle. 
Blackbody  calibration  techniques  are  used  to  allow  conversion  of  this 
radiative  signal  into  spectral  radiance  data,  which  may  then  be 
analyzed  to  provide  spectral  emissivity.  The  original  operative  con- 
cept of  the  apparatus  necessary  to  make  these  spectral  emissivity 
measurements  is  shown  in  Figure  IV-1. 

Calculations  performed  to  predict  the  temperature  achieved  by 

laser  heating  of  micron-size  AI2O2  particles  (discussed  in  detail  in 

Section  IV-B.3.C)  using  the  representative  values  of  n and  k from  the 

POPM  literature  survey  report^^^^'  indicate  that  laser  energy 

4 2 

density  fluxes  near  10  watts/cm  are  required  to  heat  4nm  to  8u,m 

diameter  particles  to  1000  K.  Figure  IV-2  presents  this  computed 

equilibrium  temperature  as  a function  of  laser  energy  flux  density. 

The  time  required  to  attain  equilibrium  is  given  in  Figure  IV-3  as  a 

function  of  laser  energy  flux,  and  is  seen  to  be  on  the  order  of  a few 

tens  of  seconds  for  4pm  to  8pm  particles  of  A1-0-.  Maximum  flux 

5 2 ^ ^ 

densities  of  1.6  X 10  watts/cm  in  a 48pm  spot  size  diameter  have 
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BACKGROUND 


Figure  IV-1.  Particle  emlsslvlty  experimental 


LASER  ENERGY  FLUX  (W/cm  ) 


been  observed  using  a single  4-watt  TEM^^  Nd:YAG  lasers  and  focus- 
ing optics  (Section  IV-B.3.b).  This  indicates  the  feasibility  of  heating 
Al^Os  particles  to  temperatures  in  excess  of  1000 °K.  Based  on 
these  results  we  expect  that  MgO,  ZrO,  and  carbdn  can  also  be  heated 

o 

to  at  least  1000  K,  since  these  materials  are  less  transparent,  i.  e.  , 
more  absorbing,  to  the  1.06|im  laser  radiation  than  is  AI2O2. 

Signal  to  noise  calculations  assuming  a collection  solid  angle  of 
rr/9,  a . 5g,m  diameter  particle,  and  an  interferometer  efficiency  of 
0.4  (Section  IV-B.4),  along  with  reasonable  detector  electronics  and 
background  assumptions,  predict  for  a 2sec  observation  time  that  17 
interferograms  are  required  to  obtain  a S/N  = 100  from  observations 

O 

of  a 900  K particle.  Thus  we  conclude  that  the  actual  detection  and 
analysis  of  radiation  from  a small  1000°K  particle  should  be  possible. 

The  levitation  of  particles  in  the  presence  of  forces  produced 

by  the  intense  radiation  from  the  heating  laser  beams  has  required 

considerable  effort.  The  commercially  purchased  Science  Spectrum 

particle  levitator  which  we  expected  to  use  proved  to  have  unreliable 

electronics,  and  a sample  cell  which  had  to  be  modified  for  the  present 

experiments.  With  the  modified  cell  and  electronics,  particles  of  all 

four  materials  of  interest  were  successfully  levitated  when  introduced 

into  the  sample  cell  by  means  of  an  aspirator  through  the  electrical 

discharge  from  a Tesla  coil.  We  found,  however,  that  a laser  energy 

2 

flux  of  50  watts/cm  , corresponding  to  the  unfocused,  polarized  out- 
put of  one  Nd:YAG  laser  exerted  sufficiently  strong  forces  (radiation 
and/or  photophoresis)  upon  the  particle  to  drive  it  from  the  potential 
well  of  the  levitator.  The  geometrical  arrangement  of  three  symmetri- 
cally placed  beams  proved  to  be  ineffective  in  preventing  particle  loss, 
although  with  modified  levitation  electronics  we  were  able  to  achieve 
retention  times  on  the  order  of  lOsec  at  full  power.  Retention  of  the 
irradiated  particle  at  the  levitation  cell  equilibrium  position  has  only 
recently  been  achieved  by  replacing  the  three  symmetrically  oriented 
laser  beams  system  with  a two-laser  opposed  beam  system  in  which 
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the  focal  points  of  the  beams,  and  the  laser  powers,  are  carefully 
adjusted. 


B.  EMISSIVITY  EXPERIMENT  DESCRIPTION 

1.  Introduction 

The  experimental  apparatus  shown  in  Figure  IV- 1 modified  to 
the  two  opposing  laser  configuration  is  shown  schematically  in 
Figure  IV-4.  Up  to  the  conclusion  of  the  work  reported  here  no 
spectral  measurements  have  been  completed,  so  the  external  infra- 
red optics  consist  of  a collection  mirror,  chopper,  and  detector 
rather  than  the  combination  of  blackbodies  and  interferometers 
planned  in  Figure  IV-4.  Two  opposed  4-watt  TEM^^  NdrYAG  lasers 
are  focused  to  a point  at  the  particle  equilibrium  position  in  the  levi- 
tator,  for  the  purpose  of  particle  heating.  Through  this  point  also 
passes  the  beam  of  an  Ar^  laser,  used  for  illuminating  the  particle 
for  position  sensing.  Figures  IV-5  and  IV -6  are  photographs  of  the 
levitation  cell  and  surrounding  apparatus.  Not  shown  are  the  Nd:YAG 
lasers  and  the  Ar^  laser.  The  rest  of  this  section  contains  details  of 
the  subsystems  which  are  critical  parts  of  the  particle  emissivity 
experiment,  as  well  as  analysis  methods  for  the  reduction  of  infrared 
emission  observations  to  obtain  emissivity. 

2.  Levitation  Apparatus 
a.  Cell  Design 

The  configuration  of  the  Science  Spectrum  Differential  II 
levitation  cell  supplied  by  the  manufacturer,  shown  in  Figure  IV- 7, 
has  been  used  except  for  modifications  to  the  transparent  window 
section.  This  portion  of  the  cell  has  been  modified  to  accept  thin 
quartz  windows  to  pass  the  Nd:YAG  heating  laser  beams.  An  IR  trans- 
parent window  for  viewing  the  thermal  radiation  from  the  particle  was 
also  added.  It  was  found  necessary  to  use  5- 15cm  long  6mm  diameter 
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Figure  IV-6.  Close-up  of  the  levitation  cell,  focusing  lenses  and  telescope.  The  chopper  and 

collection  mirror  arc  on  the  right  far  side  of  the  levitation  cell.  The  beam  expanders 
for  the  Nd:YAG  lasers  are  shown  on  the  left. 
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Figure  IV-7.  Scattering  cell. 


tubes,  painted  black  inside,  to  space  the  Nd;YAG  and  argon  ion  laser 
beam  windows  away  from  the  cell.  This  greatly  reduced  the  levels 
of  scattered  light  in  the  cell,  and  thus  increased  the  sensitivity  of  the 
optical  servo  system. 

Particles  v/ere  charged  in  an  RF  discharge  located  in  a glass  tube 
through  which  sample  particles  were  blown  by  a bulb.  These  particles 
passed  through  a length  of  rubber  tubing,  into  the  setting  chamber,  and 
on  into  the  cell  with  sufficient  charge  for  levitation.  The  exhaust  con- 
nector was  sealed  closed  after  particle  injection  to  prevent  variations 
in  outside  air  pressure  from  displacing  a levitated  particle. 

b.  Electronics 

Several  months  of  operation  using  the  Science  Spectrum 
Differential  II  levitation  electronics  demonstrated  that  these  electron- 
ics were  unreliable  and  difficult  to  repair.  Interference  problems 
from  the  RF  discharge  unit  were  insoluble,  and  required  that  the 
levitator  be  turned  off  while  charged  particles  were  injected.  In 
addition,  the  maximum  voltage  output  for  the  unit  was  ^ 300  volts, 
which  was  not  high  enough  to  generate  a potential  well  within  the  levi- 
tation cell  of  sufficient  strength  to  hold  particles  in  the  presence  of 
radiation  produced  forces.  A new  set  of  levitation  electronics  was 
therefore  built.  The  block  diagram  of  this  unit  is  given  in  Figure  IV-8. 
The  image  of  the  levitated  particle  illuminated  by  the  argon  ion  laser 
is  split  horizontally  and  projected  onto  two  4818  photomultipliers. 

We  retained  the  original  optics  designed  for  this  purpose,  but  added 
a 6lcm  long  1:1  telescope  to  space  the  electro-optical  servo  box 
from  the  levitation  cell.  In  addition,  the  two  photomultiplier  high 
voltage  circuits  were  separated,  so  fine  adjustments  in  the  photo- 
multiplier voltages  could  be  made  to  compensate  for  slight  differences 
in  tube  gains.  Type  4818  tubes  were  employed  instead  of  the  original 
1P21  tubes  since  these  have  a high  gain  with  very  low  background 
noise  level. 
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Figure  IV-8.  Block  diagram  of  levitation  electronics. 


r 
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The  two  signals  from  the  photomultipliers  are  processed  by  a 

log -difference  circuit,  and  this  output  signal  ultimately  determines 

the  voltages  and  in  the  cell.  The  top  plate  of  the  cell  is 

ground,  and  V . is  2 x V , : the  actual  signs  of  V . and  V , • 

pin  plate  pm  plate 

depend  upon  the  charge  of  the  levitated  particle.  Two  Ortec  Model  456 

power  supplies,  slightly  modified  for  msec  voltage  tracking  character-  t 

istics,  are  used  as  high  voltage  supplies.  The  highest  voltages 

attainable  are  on  the  order  of  3kV  but  such  high  voltages  cause  arcing  I 

within  the  cell,  so  actual  operation  is  usually  in  the  600  to  800  volt  ' 

range  for  V . . I | 

pm  ' 

A future  planned  modification  for  this  servo  system  involves 

the  replacement  of  the  two  photomultipliers  by  a single  split  field 

solid  state  detector,  UDT  PIN  SPOT/8D.  This  will  eliminate  two 

high  voltage  power  supplies,  a bulky  optical  system,  and  an  expensive 

argon  ion  laser.  The  laser  will  be  replaced  by  a low-cost  5mW  HeNe 

unit. 


c.  Operational  Details 

Successful  operation  of  the  levitator  is  a skill  which 
requires  several  hours  of  practice  to  acquire;  however,  particles  can 
be  captured  and  levitated  in  a matter  of  minutes  by  an  experienced 
operator.  The  general  procedure  for  levitating  particles  is  as  follows: 

(1)  Place  the  levitator  electronics  in  the  manual  mode.  This 
disconnects  the  optical  servo  and  allows  the  operator  to 
inanually  adjust  voltages  to  hold  a particle  in  the  levitator 
during  the  period  of  time  when  other  particles  are  settling 
out  of  the  optical  system  field  of  view. 

(2)  Set  the  system  gain  to  a level  which  results  in  only  100  to 
200  volts  maximum  for  V . . Higher  gains  make  manual 
operation  very  difficult  and  result  in  the  capture  of  parti- 
cles with  a smaller  charge. 
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(3)  Turn  on  the  RF  discharge  and  blow  a cloud  of  charged 
particles  Into  the  cell. 

(4)  Capture  a particle  of  the  appropriate  charge  and  hold  it 
in  position  by  manually  adjusting  the  voltages  until  all 
other  particles  have  settled  out  of  the  cell. 

(5)  Turn  the  levitator  electronics  to  the  servo  mode;  the 
particle  should  be  held  automatically  at  the  equilibrium 
position. 

(6)  Increase  the  gain  to  a value  which  results  in  V . between 

pin 

600  and  800  volts  maximum;  this  should  strongly  hold  the 
particle  and  greatly  reduce  any  oscillations. 

3.  Particle  Irradiation  and  Heating 

The  optical  configuration  required  for  successful  irradiation  of 
a levitated  particle  has  been  given  in  Figure  IV-4.  Two  focused, 
opposed  Nd:YAG  lasers  are  used,  together  with  coaxial  HeNe  align- 
ment lasers.  The  laser  output  is  first  expanded  using  a Holobeam  type 
1006190  upcollimator,  then  focused  to  a point  approximately  48pm  in 
diameter  by  a 2.  5cm  diameter,  20cm  focal  length  lens  made  from  fused 
silica  and  antireflection  coated  for  use  at  1.06p,m.  This  focusing  lens 
has  a depth  of  field  on  the  order  of  5mm.  The  lasers  used  for  heating 
are  two  Holobeam  255  Nd;YAG  CW  lasers  used  in  the  TEMqq  mode. 
Typical  single  mode  output  is  ~4  watts  CW.  The  laser  energy  flux 
measured  at  the  focal  point  in  the  cell  is  ~ 1.  6 X 10^  watts/cm^.  This 
is  sufficient  to  heat  particles  in  the  size  range  of  interest  to  temper- 
atures in  excess  of  1000  K,  as  will  be  shown  shortly.  For  convenient 
alignment  purposes  a Spectra-Physics  Model  156  Imwatt  HeNe  laser 
is  mounted  with  the  Nd;YAG  lasers  and  aligned  to  produce  a .63p,m 
beam  colinear  with  the  l,06|jm  Nd:YAG  beam. 
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a.  Alignment 

Alignment  of  the  laser  beams  is  a critical  problem  in  this 
experiment  because  the  Nd:YAG  and  levltator  beams  must  intersect 
within  a volume  less  than  ZOpm  in  diameter.  We  have  demonstrated 
that  the  use  of  Oriel  Model  1460  high  precision  mounts  for  beam 
adjustment  mirrors  allows  micron-scale  adjustments  within  the  levi- 
tation cell.  Initial  beam  adjustments  are  performed  using  the  HeNe 
alignment  laser.  Final  adjustments  are  then  made  with  the  Nd:YAG 
scattered  signal  from  the  particle,  using  an  IR  viewer. 

Experiments  have  shown  that  a single  focused  NdrYAG  laser 
irradiating  a particle  in  the  levitator  rapidly  moves  the  particle  away 
from  the  equilibrium  position.  The  original  concept  of  balancing 
laser-induced  forces  on  the  particle  by  means  of  three  symmetrically 
placed  beams  has  proven  impossible  to  implement.  Particle  residence 
times  in  excess  of  lOsec  in  the  cell  were  achieved  with  high  cell 
voltages  and  some  minor  changes  in  cell  electrode  configuration,  but 
the  particle  was  ultimately  lost. 

The  configuration  of  two  opposed  lasers  has  proven  to  be  success- 
ful in  irradiating  the  particle  without  moving  it.  Adjustment  of  the 
laser  power  levels  to  the  point  where  the  forces  effectively  cancel  has 
proven  to  be  quite  difficult,  and  typically  occupies  approximately  30 
minutes.  With  the  irradiated  particle  in  a stable  position,  observa- 
tions on  the  order  of  hours  are  possible. 

The  focal  points  of  the  two  lasers  must  be  carefully  placed  to 
entrap  the  levitated  particle,  as  shown  in  Figure  IV-9.  The  radiative 
forces  are  proportional  to  the  flux  density  in  laser  beams.  Thus,  for 
example,  if  the  particle  moves  to  the  left  F2  increases  and  F^  decrea- 
ses, tending  to  return  the  particle  to  the  neutral  point.  In  fact,  it  is 
unimportant  which  direction  the  laser  induced  forces  act  on  the  parti- 
cle since  an  adjustment  of  the  focusing  optics  to  interchange  the  posi- 
tions of  the  focal  points  of  the  two  lasers  will  compensate  if  the  forces 
act  in  the  opposite  direction.  This  opposed  laser  balancing  is  required 
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opposed  laser  configuration. 


because  we  observe  that  the  radiation  forces  are  strong  compared  to 
the  potential  well  of  the  levitator.  These  observations  are  consistent 
with  the  order  of  magnitude  calculations  for  photon  momentum  trans- 
fer to  a particle  through  either  absorption  or  reflection,  and  is  illus- 


trated in  Figure  IV- 10. 


b.  Laser  Flux  Density  Measurements 

The  success  of  the  particle  emissivity  experiment  in  the 
POPM  program  depends  critically  on  being  able  to  heat  the  suspended 

particle  with  a high  flux  density  laser  beam.  For  micron-size  ALO„ 

5 2 ^ ^ 

particles  flux  densities  of  ~ 10  watts/cm  on  the  particle  are  needed 

to  achieve  temperatures  near  the  melting  point.  Given  a nominal 

five  watt  laser,  this  requires  that  the  spot  size,  at  the  particle,  be 

of  the  order  of  tens  of  microns  in  order  to  achieve  the  required  flux 

density. 

An  experiment  was  performed  to  measure  the  maximum  flux 
density  obtainable  from  a nominal  five  watt,  single  mode,  Nd:YAG 
laser.  A description  of  this  exneriment  and  the  results  are  contained 
in  this  section. 

The  flux  density  contour  for  a laser  operated  in  the  TEMqq^ 
mode  is  Gaussian,  and  may  be  written 


: i 


1 


1 


, 2 2,,  2 2,2 
F(r)  = F„e-<"  ^ = F^e”- 

2 2 1/2 

where  F(r)  is  the  flux  density  at  a distance  r = (x  + y ) away  from  j 

the  beam  axis,  F„  = F(0),  and  a is  a measure  of  how  narrow  the  beam 

^ - 1 i 

is.  For  instance,  when  r = a,  F(a)/FQ  = e , the  so-called  "one  over  | 

e point".  A plot  of  F(r)  as  a function  of  x and  y (or  r)  is  shown  in  j 

Figure  IV-11.  This  is  the  Gaussian  or  "bell  shaped"  surface,  where  | 

X and  y are  the  orthogonal  spatial  coordinates  perpendicular  to  the  i 

laser  beam  axis  and  the  z coordinate  gives  F(r),  the  flux  density,  at  ^ 
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Figure  IV-10.  Computed  acceleration  of  a 5p,m  carbon  particle  due  to  radiation  pressure  and 
neglecting  all  drag  forces. 


any  distance  r from  the  laser  beam  axis.  We  now  calculate  the  power 
that  would  be  passed  by  a circular  aperture  of  radius  r centered  on 

3. 

the  laser  beam  axis* 


P(r^)  = \ F(r)dA  = \ ^ \ F^e'""  rdrdV 

^ -'a  >^0  -0  ° 

2 2 

The  integration  over  <P  is  trivial,  and  multiplying  by  (-2a  /-2a  ) we 
have 


P(r,)  = - 

and  since  the  term  in  parenthesis  is  the  differential  of  the  argument 
of  the  exponential,  we  have 


P(ra) 


Letting  r — ►“  gives  the  total  power  in  the  beeim,  P(total),  that  is 
a 

P(total)  = to^Fq 


Therefore,  the  ratio  of  the  power  passed  with  the  circular 
aperture  to  the  total  power  is. 


P(total) 


(IV- 1 


Thus  a measure  of  this  ratio  with  a known  pinhole  radius  determines 
a,  or  how  narrow  the  beam  is. 
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A schematic  of  the  laser  flux  density  experiment  is  shown  in 
Figure  IV-12.  The  laser  beam  is  expanded  and  goes  into  a beam 
alignment  device  which  allows  the  beam  to  be  translated  in  any  direc- 
tion perpendicular  to  the  axis  of  the  beam.  The  beam  is  focused  by 
a lens  and  somewhat  after  the  focus  point  is  located  a detector  which 
intercepts  the  entire  beam.  The  beam  profile  is  probed  by  means  of 
a pinhole  aperture  which  can  be  moved  along  the  axis  of  the  beam. 

The  beam  diameter,  o,  will  be  a function  of  where  along  the 
beam  the  pinhole  is  located,  and  by  definition  has  a minimum  value  at 
the  focus.  Flux  density  will  be  a maximum  where  a is  a minimum. 

The  experimental  data  were  obtained  as  follows.  With  the  pin- 
hole aperture  removed,  the  laser  was  adjusted  to  give  three  watts  at 
the  detector  (detector  element  area  large  with  respect  to  beam  dia- 
meter). The  pinhole  aperture  was  mounted  on  a micrometer  stage 
such  that  its  12mm  travel  was  approximately  centered  on  the  focal 
point.  Power  passed  by  the  pinhole  aperture  was  measured  as  a 
function  of  distance  (d  in  Figure  IV-12)  in  steps  of  0.25mm.  The  pin- 
hole aperture  center  could  not  be  made  to  center  on  the  laser  beam 
axis  exactly  (crude  analysis  indicates  it  moved  about  20^  over  a 
travel  of  3.  75mm).  However,  the  beam  alignment  device  provided 
the  mechanism  whereby  the  beam  could  be  translated  so  that  it  could 
be  cantered  on  the  pinhole.  A typical  set  of  data  is  given  in  Table  IV -1. 
The  first  column  gives  the  micrometer  reading,  and  the  second  column 
the  number  of  watts  passed  by  the  99tim  diameter  pinhole  aperture. 

The  third  column  is  the  ratio  P(r  )/P(total),  and  the  fourth  column  is 
the  value  of  a calculated  from  equation  IV- 1.  The  data  are  plotted  in 
Figure  IV-13.  The  error  bars  were  determined  assuming  an  error  of 
i 0.  05  watts  and  are  representative  of  the  noise  encountered  in  the 
measurements.  A sm.ooth  curve  was  drawn,  and  although  it  is  not  a 
fit  in  the  least  squares  sense,  it  is  estimated  from  the  minimum  in 
the  curve  and  the  size  of  the  error  bars  that. 


J 
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Figure  IV- 12.  Schematic  of  focused  heating  laser  flux  measurement  apparatus. 


Table  IV- 1.  Pinhole  probe  data  of  laser  Intensity, 


Power 

P(ri)/ 

d(mm) 

(watts) 

P(total) 

o(|j,rr 

6.  50 

2.45 

0.  817 

38.  C 

6.  75 

2.  65 

0.  883 

33.  i 

7.  00 

2.  82 

0.  940 

29.  E 

7.  25 

2.  85 

0.  950 

28.  ( 

7.  50 

2.  87 

0.  957 

27.  S 

7.  75 

2.  90 

0.  967 

26,  i 

8.  00 

2.  92 

0.  973 

26.  ] 

8.  25 

2.  97 

0.  990 

23.  1 

8.  50 

2.  95 

0.  983 

24.  E 

8.  75 

2.  90 

0.  967 

26.  8 

9.  00 

2.  85 

0.  950 

28.  8 

9.25 

2.75 

0.  917 

31.4 

9.  50 

2.  70 

0.  900 

32.  ( 

9.  75 

2.45 

0.  817 

38.  C 

P(total  = 3.  00  watts 
r = 49.  S^m 


IV- 13.  Heating  laser  beam  diameter  a versus  d from  pinhole  data 


l»S£SSS 


When  r = a. 


a(xnin)  = 24  ± 2|j,m. 


- 1 -e-‘~  0.6321, 


which  says  63.  21  percent  of  the  power  in  the  beam  is  confined  in  a 
circle  whose  radius  is  a.  The  average  flux  density  is  then, 


p,  > _ (0.  6321)  P(total)  watts 

' ave  2 2 


When  the  laser  is  optimally  tuned  and  the  excitation  lamp  is  operated 
near  maximum  output,  P(total)  =>■  4.  5 watts  and  using  a = 24p,m,  we 
obtain 


F(a)  = 1.  57  X 10^  watts/cm^ 
' avg 


The  choice  of  the  "one  over  e point"  as  the  limiting  radius  is  some- 
what arbitrary.  What  it  means,  however,  is  that  if  we  could  center 

a 48gm  diameter  particle  in  the  laser  beam,  it  would  experience  an 

5 2 

average  flux  density  of  ~ 1. 6 X 10  watts/cm  . For  a smaller  parti- 
cle the  average  flux  would  be  higher.  For  instance,  for  a 24p,m  dia- 
meter particle  (i.  e. , a/2)  we  have 


F(a/2)  = (0.22121(4.5)  w^  ^ 2.20  x 10^ 

tt(1.2x10  ) cm  cm'^ 


It  will  be  shown  in  the  following  section  that  it  should  be  possible  to 
melt  Al^Oj  micron-sized  particles  with  flux  densities  in  this  range. 
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Figure  IV- 13  also  provides  an  estimate  of  the  depth  of  field  we 
have  at  our  disposal.  If  we  allow  a to  vary  by  ~ 10  percent  (2.  5p,m) 
the  allowable  variation  in  d is  ~ 1.3mm  (see  Figure  IV-12).  That  is, 
if  we  locate  the  particle  to  within  ± 0.  65mm  of  the  true  a(min)  we  will 
have  a maximum  degradation  of  ~ 20  percent  in  flux  density. 

c.  Laser  Heating  of  Particles 

This  subsection  describes  the  results  of  computer  calcu- 
lations of  the  laser  particle  heating  scheme  for  aluminum  oxide. 
Ignoring  convection  and  conduction  the  basic  heat  transfer  equation 
applicable  to  laser  irradiation  heating  is: 


Hct  = |nr^ 


^ -5  -5 

dT  r ^1^  p ^1^ 

"S  dT  + i iSScTxrrrr  '»■  - S?(c,/tTJ-T 


CjX 

2 


where  H is  the  incident  laser  power,  a the  absorption  efficiency  factor, 
r the  particle  radius,  p the  particle  density,  the  heat  capacity  of 
the  particle,  T the  temperature,  Tq  the  temperature  of  the  surround- 
ings, t the  time,  e the  emissivity,  the  first  radiation  constant, 
the  second  radiation  constant,  and  X the  wavelength  of  the  radiation 
emitted.  The  term  on  the  left  describes  the  steady-state  energy  input 
to  the  particle  from  the  laser.  The  first  term  on  the  right  side  of  the 
equation  describes  the  heat  stored  within  the  particle.  The  second 
term  describes  the  energy  emitted  by  the  particle  as  a blackbody. 

The  last  term  describes  the  energy  absorbed  by  the  particle  from  the 
surroundings. 

A program  was  written  which  computed  the  particle  temperature 
as  a function  of  time.  The  calculations  were  carried  out  long  enough 
in  time  for  equilibrium  to  occur.  Since  equation  IV-2  is  a differential 
equation  in  time,  the  methods  of  numerical  analysis  had  to  be  used  to 
perform  the  integration.  The  Runge-Kutta  method  was  used  to  perform 
the  integration. 


-180- 


The  parameters  of  equation  IV-2  were  evaluated  from  the  best 
available  data  for  aluminum  oxide.  Three  values  of  the  particle 
radius  were  used  for  the  calculations  — 0.  5,  4.  0 and  8.  0|ji,m.  The 
density  of  was  taken  to  be  3.  965gms /cm^.  The  heat 

capacity  of  aluminum  oxide  is  a function  of  temperature  as  is  shown  in 
Figure  IV-14.^^®^'  These  data  were  tabulated  at  closely  spaced 

temperature  intervals  and  interpolated  linearly  by  the  computer.  The 
emissivity  of  aluminum  oxide  was  taken  to  be  0.  85  in  the  wavelength 
interval  of  4.  5 to  11.  5p,m  and  0.  0 outside  this  range.  This  is  a rea- 
sonable approximation  of  the  data  of  Mergerian^^®^'  which  is 

shown  in  Figure  IV-15.  The  blackbody  integrals  were  evaluated  over 
the  corresponding  wavelength  interval  — 4.  5 to  11.5^,m.  The  effic- 
iency factor  for  absorption,  a,  was  calculated  from  Mie  theory  for 
X = 1.06p,m  and  the  three  particle  radii  given  above  from  absorption 
coefficient  data.  The  absorption  coefficient  of  aluminum  oxide  is 
also  a function  of  temperature.  Unfortunately,  there  does  not  exist 
data  by  the  same  experimenter  and  apparatus  over  the  whole  temper- 
ature range  of  interest.  Mergerian^^^^’  covers  the  range  from 

373  to  1273  K.  Gryvnak  and  Burch^^^^*  cover  the  range  from 

O 

1473  to  1973  K.  Unfortunately,  the  two  sets  of  data  do  not  join 
smoothly.  The  calculated  efficiency  factors  for  absorption  are  shown 
in  Figure  IV- 16.  The  discontinuity  around  1400 °K  is  obvious.  The 
data  were  linearly  interpolated  by  the  computer.  For  some  high 
temperature  calculations  an  artificial  value  of  a at  1473  °K  was  inserted. 
This  artificial  value  is  a simple  interpolation  of  the  data  at  1273  and 

O 

1773  K and  are  shown  as  solid  symbols  on  Figure  IV- 16.  The  laser 

5 2 

power  was  varied  between  40  and  4X10  watts/cm  . The  laser 
intensity  measurements  described  in  the  previous  subsection  show 
that  the  Nd;YAG  lasers  used  in  the  particle  experiment  are  capable  of 
fluxes  as  high  as  2 X 10  watts /cm  . 
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HEAT  CAPACITY,  C (Joules /gm  KI 


Figure  IV-14.  Heat  capacity  of  aluminum  oxide,  AI2O 
(from  Lange,  Ref.  IV-<i). 


WAVELENGTH  IN  MICRONS 
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The  results  of  the  calculation  are  given  in  Table  IV-2  and  in 
Figures  IV-2  and  IV-3.  Figure  IV-2  shows  the  equilibrium  tempera- 
ture attained  as  a function  of  laser  power.  Figure  IV-3  shows  the 
time  required  to  achieve  this  equilibrium  temperature  as  a function 
of  laser  power.  The  open  symbols  are  for  the  unsmoothed  a data  and 
the  solid  symbols  for  the  artificially  smoothed  a data. 

The  results  shown  in  Table  IV-2  show  an  interesting  feature  of 
the  time  history  of  the  heating  of  a particle.  That  is  the  asymptotic 
approach  to  equilibrium.  The  length  of  time  to  achieve  the  last  one 
percent  of  the  equilibrium  temperature,  particularly  at  low  laser 
power  levels,  can  be  significant.  For  example,  16pm  diameter 
particles  with  40  watts/cm  incident  laser  power  attain  99  percent  of 
the  equilibrium  temperature  in  less  than  one  minute  but  require 
almost  three  more  minutes  to  achieve  the  last  one  percent.  Great 
care  must  be  exercised  when  conducting  the  particle  emissivity  exper- 
iments to  allow  sufficient  time  for  equilibrium  to  be  attained.  Figures 
IV- 17  and  IV- 18  show  the  time  history  of  the  heating  of  various  dia- 
meter particles  for  two  different  laser  power  levels.  In  general,  the 
small  particles  equilibrate  much  more  rapidly  than  the  larger  parti- 
cles. Also,  the  higher  the  laser  power,  the  faster  equilibrium  is 
attained. 

The  calculations  indicate  a potential  problem  with  aluminum 
oxide.  This  problem  is  graphically  illustrated  in  Figure  IV-2  with 
the  step  function  behavior  of  the  equilibrium  temperature.  The  prob- 
lem is  thermal  runaway  caused  by  the  precipitous  increase  in  the 
efficiency  factor  for  absorption  of  energy  above  1500  K.  Figure  IV- 19 
is  an  additional  illustration  of  this  problem.  This  figure  shows  the 

time  history  of  the  heating  of  a 16pm  diameter  particle  with  an  inci- 

4 2 

dent  laser  power  level  of  6 X 10  watts/cm  . The  calculations  were 
made  with  the  smoothed  data  for  the  efficiency  factor  for  absorption. 

O 

The  temperature  of  the  particle  is  well  behaved  up  to  1600-1700  K at 
which  point  thermal  runaway  occurs.  If  the  absorption  coefficient  of 
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Al^Os  Is  truly  as  Indicated  in  Figure  IV- 16,  there  will  be  great  diffi- 
culty in  achieving  stable  temperatures  above  1600°K  since  the  particle 
temperature  appears  to  be  extremely  sensitive  to  small  variations  in 
laser  power.  However,  the  absorption  coefficient  is  not  known  well 
as  a function  of  temperature  and  the  runaway  problem  may  not  be  as 
severe  as  these  calculations  would  indicate. 

These  calculations , together  with  the  measured  laser  intensities 
reported  in  the  previous  subsection  indicate  that  heating  particles  to 

O 

at  least  1000  K should  not  be  difficult  in  principle  if  the  losses  due  to 
convection  and  conduction  are  not  large.  If  these  additional  losses 
are  large  they  can  be  decreased  by  evacuating  the  cell,  which  will 
require  changes  in  the  damping  of  the  particle  levitation  servo  system. 

4.  Infrared  Radiation  Observation 

In  this  subsection  signal-to-noise  calculations,  particle  temper- 
ature measurement  calculations,  and  optical  configurations  will  be 
described.  We  had  the  opportunity  to  make  only  two  attempts  to 
observe  radiation  in  the  few  days  between  the  success  of  the  stable 
two-laser  configuration  and  the  relocation  of  the  POPM  laboratory  to 
a new  building.  Both  of  these  rather  quickly  assembled  attempts 
were  unsuccessful,  we  believe  due  to  the  fact  that  the  stringent  posi- 
tioning requirements  of  the  infrared  detector  relative  to  the  image  of 
the  particle  were  not  met. 

The  present  optical  system  consists  of  several  elements.  Rad- 
iation from  the  particle  passes  through  a . 75  x 1.  5cm  window  in  the 
cell  wall,  and  is  collected  by  a 7.  95cm  diameter,  23  degree  off-axis 
paraboloidal  mirror,  which  focuses  the  radiation  at  a distance  of  - Im 
from  the  cell.  Placed  between  the  cell  and  the  mirror  is  a 790Hz  cir- 
cular chopper.  The  radiation  is  detected  by  a liquid  helium-cooled 
copper  doped  germanium  detector,  Raytheon  QKNl  546.  The  detector 
characteristics  are:  KRS-5  entrance  window,  7.85  x lO'^cm^  area, 
five  degrees  field-of- view,  andN.E.P.  of  6.  8 x 10 ^ WHz~  ^ The 
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procedure  for  optical  alLgriment  will  be  to  use  the  detector  without  a 
blocking  filter  to  detect  the  scattered  l.Obyjn  Nd;YAG  laser  radiation 
from  the  particle.  When  this  is  accomplished,  a broadband  filter  will 
be  placed  in  front  of  the  detector,  and  the  radiation  from  the  heated 
particle  will  be  examined. 

The  final  optical  system  is  shown  in  Figure  IV-4.  The  same 

collection  mirror  will  be  employed,  the  converging  particle  radiation 

will  pass  through  a Digilab  fast  scanning  Mlchelson  interferometer 

and  be  chopped  at  ~ 9kHz  using  a high  speed  rotating  chopper.  Laser 

Precision  CTX  534.  The  detector  to  be  used  will  be  a Raytheon 

-4  2 

copper-doped  germanium  wi •'n  area  of  5 x 10  cm  and  2.5  degree 
field  of  view.  The  subseq'  nal  processing  will  be  similar  to 

that  already  described  for  th.  oulk  experiment. 

a.  Signal-to-Noise  Calculations 

The  hot  particles  are  small  in  size  and  therefore  provide 
only  a small  infrared  signal.  In  this  subsection  we  show  by  calcula- 
tion that  adequate  signal  levels  should  be  available  to  make  spectral 
measurements  of  the  particle  radiance  with  the  multiplex  advantage 
of  Fourier  transform  spectroscopy.  To  carry  out  this  calculation 
we  must  begin  by  summarizing  the  performance  characteristics  of 
our  interferometer. 

The  Aerospace  Michelson  inlc  rfc  romcler  can  operate  in  the 
wav^elengtli  region  ~2.8  to  25um  (v  ~ 3 600  to  400cm  S.  The  stationary 
mirror  of  the  interferometer  has  been  displaced  0.  25cm)  so  that  two 
sided  interfercjgrams  (-  0.25cm  ^ x 0.25cm)  can  be  obtained  which 
can  be  transformed  into  spectra  witliout  the  need  for  phase  correction. 
The  corresponding  resolution  is  4cm  \ Howexer,  with  little  diffi- 
culty, the  stationary  mirror  could  be  further  displaced,  so  that  a 
spectral  resolution  of  ~ 1cm  could  he  achiexed  if  desired. 


■ds. 


The  scan  velocity  of  the  moving  mirror  Is  ^ 0.  5 optical  cm/ sec 

which  for  V = 400cm" ^ yields  f^  = 200Hz.  (That  Is  a spectral  "line" 

- 1 ^ 

at  V = 400cm  would  have  as  Its  Interferogram  a cosine  wave  of 
frequency  200Hz.  ) Likewise,  for  v = 3600cm  \ f^^  = 1800Hz  for  a 
nominal  bandwidth  of  ~ l600Hz(=  Af).  This  bandwidth  will  have  to  be 
Increased  slightly  because  of  the  practical  limitations  of  the  Inter- 
ferometer. 

Sampling  tlieory  dictates  that  the  Interferogram  must  be  sam- 
pled at  least  as  frequently  (In  optical  path  difference)  as 


Ax  = 


where  v Is  the  largest  wavenumber,  with  non-zero  Intensity,  In  the 
^ - 1 

spectrum.  In  our  case  ~ 3600cm  , due  to  optical  components  In 

the  Interferometer,  so  Ax(max)  ~ 1.3889urn.  However,  Ax  for  our 
Interferometer  Is  measured  In  Integral  powers  of  two  times  half  the 
wavelength  of  the  HeNe  laser  line  (X  = 0.6328um),  that  Is, 


Ax  = j,  X,  2X,  4X,  8X,  etc.  For  our  case. 


/ix(optlmum)  _ 1.3889fii^'>  _ ^ j 
X/2  " 0.3l64pim  " * 


Therefore,  Ax(practlcal)  = 1.2656  x 10  cm.  The  number  of 
points  per  interferogram  n,  is  then. 


optical  path  traversed  _ , , ,^3  . 

" = -rxlpra^tlEII) P'®- 


The  fast  Fourier  transform  program  efficiently  transforms 

m 3 

intcrferograms  with  n = 2 . For  m = 12,  n = 4.096  « 10  pts.  Work- 

ing backwards  and  remembering  that  two-sided  intcrferograms  are 
l)i‘lng  measured. 
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L = ^ = (2048)  (1.2656  x 10  ) = 0.2592cm, 

which  yields  a resolution  of  ~ 3.  86cm  Also,  Ax(practlcal)  = 

-4  - 1 

1.2656  X 10  cm,  and  so  3950cm  , which  for  our  moving  mir- 

ror velocity  of  0.  5cm/sec  yields  f„  = 1975Hz  and  Af  = 1775Hz.  The 
region  from  ~ 3600cm”  to  3950cm"^  will  be  "empty"  (i.  e.  , no 
spectra)  because  of  the  natural  cut-off  at  v 3600cm”^.  Although 
there  will  be  no  useful  spectroscopic  data  in  this  region  it  can  be 

used  to  estimate  the  noise  in  the  "occupied"  spectral  region  as  has 

[Ref  IV-51 

been  pointed  out  by  Connes.*^  ' ” ■'  This  does  not  imply  that  only 

this  method  is  used  to  measure  system  noise. 

It  is  now  possible  to  proceed  with  the  signal  to  noise  computa- 
tion. Assume  a background  limited  detector  on  the  interferometer 
with  D='=  = 5 X lO^^cm  Hz^^^  W a n/ 9 field  of  view  and  an  area 
A = 10  Then 

NEP  = = 2.66  X 10”^V 

where  we  have  assumed  a Af  = 1775Hz  from  the  spectral  range  of 
interest  and  the  parameters  of  our  interferometer.  It  has  been  shown 

that  the  noise  in  the  spectrum,  N , is  related  to  the  noise  in  the  inter- 

. [Ref.  IV-6]  , ® 

ferogram  ^ by. 

Ns  = ~ 1- 06  X 10‘^^(W/cm) 

- 1 -4 

where  L = 0.  25cm  (resolution  = 4cm  ),  Ax  = 1 . 2656  X 10  cm  and 
is  the  NEP  as  calculated  above.  The  above  expression  is  appropriate 
for  a single-sided  inte rfe rogram.  For  a double-sided  interferogram 
we  pick  up  a factor  of  \/~Z,  and  furthermore  the  best  resolution 
required  is  about  three  resolution  widths  so  we  can  integrate  o\er 
these  three  elements  and  pick  up  a further  factor  of  \/~3.  Thus, 


F 
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N (~  12cm  element)  = 4.  32  X 10  watts, 
s 

To  estimate  a signal  for  comparison  with  the  noise  we  assume  a 
particle  of  radius,  r = 0.5pm  and  emissivity  e = 1.  This  corresponds 
to  about  the  smallest  carbon  particle  size  to  be  investigated  in  this 
project  and  will  give  a reasonably  conservative  estimate  of  the  expected 
signal.  Furthermore,  assume  the  solid  angle  over  which  radiation  is 
collected,  Cl  = v / 9 steradians,  and  the  transmittance  or  efficiency  of 
the  interferometer  is  f = 0.4.  Then  the  detected  particle  flux  is, 


N = Aef2fN^(T), 

P B 

where  A = nr^,  the  effective  emitting  particle  area,  and  N2(T)  is  the 
blackbody  radiance  at  temperature  T.  The  above  expression  was  eval- 
uated for  a range  of  temperatures  (400  to  900  K)  for  the  spectral  ele- 
ment 400  to  412cm'  which  again  should  lead  to  a conservative  estimate 
since  for  the  tem.peratures  considered  this  wavenumber  region  lies  on 
the  longe  r wavelength  side  of  the  blackbody  radiance  function.  With 
these  assumptions  Table  IV-3  was  constructed.  For  the  various  parti- 
cle temperatures  listed  in  the  first  column,  the  second  column  gives 
the  integrated  (v  = 400  to  412cm'^)  signal-to-noise  ratio  for  a single 
interferogram.  The  third  column  gives  the  number  of  coadded  inter- 
ferograms  required  to  increase  the  signal-to-noise  ratio  to  100  in  this 
same  spectral  element;  The  fourth  column  gives  the  tinie  required  to 
measure  this  number  of  interferograms.  The  table  indicates  that  if 
the  temperature  of  the  emitting  particle  could  be  changed  quickly,  and 
the  detector  noise  were  the  limiting  noise  source,  the  data  could  be 
gathered  in  ten  minutes  or  somewhat  longer,  the  extra  time  being  used 
to  obtain  background  interferograms.  The  S/N  and  related  numbers 
given  above  are  optimistic,  especially  for  the  lower  temperatures, 
since  no  account  has  been  given  to  the  - 300  K bath  in  which  the  emit- 
ting particle  will  probably  be  immersed.  To  suppress  unwanted  sig- 

O 

nals  due  to  radiation  from  the  300  K surroundings  in  this  apparatus. 
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Table  IV-3.  Computed  s Ignal-to-nolse  for  0.  Sym  particle  at 
various  temperatures  (v  = 400  to  41Zcm 


I 

I 

/ 

1 

j 


T 

(K) 

Signal-to- 
Noise  Ratio 
one  Lnterferogram 

n.  Number  of 
Interferograms  to 
achieve  S/N  = 100 

Time  to  Obtain 
n interferograms 
(minutes) 

400 

6 

277 

4.  6 

500 

10 

100 

1.7 

600 

13 

59 

1.0 

700 

16 

39 

0.  7 

800 

20 

25 

0.  4 

900 

24 

17 

0.  3 

I 
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a modulated  signal  is  sent  to  the  interferometer  used  to  make  spectral 
measurements.  The  chopping  frequency  must  be  significantly  greater 
than  the  electrical  frequencies  corresponding  to  the  wavelength  region 
of  interest.  For  the  scan  speeds  of  our  Michelson  interferometers 
used  to  obtain  spectral  data,  the  highest  frequency  of  interest  is 
- ZkHz.  The  chopper  operates  at  ~ 8kHz  and  thus  adequate  sampling 
is  achieved.  However,  a background  problem  still  exists. 

The  infrared  image  of  the  radiating  particle  is  focused  on  a 
chopper  blade  in  the  present  appara,tus.  The  modulated  signal  con- 
sists of  the  particle  emission,  plus' emission  from  the  particle  sur- 
roundings, minus  the  emission  of  the  chopper  blade.  Ideally,  the 
emission  for  the  particle  surroundings  and  the  chopper  blade  would  be 
identical  so  that  when  no  hot  particle  is  present  and  the  chopper  is 
running,  there  will  be  no  AC  signal.  However,  if  the  blade  and  the 
background  are  at  different  temperatures,  there  will  be  a residual 
AC  background  which  must  be  distinguished  from  the  signal  due  to 
particle  emission. 


In  order  to  estimate  the  severity  of  this  background  probleni, 
the  following  analysis  was  carried  out.  Using  the  optical  paranieters 
of  the  particle  emissivity  apparatus,  a number  of  curves  were  gen- 
erated for  frequencies  between  500  and  2500cm  ^ (20  to  4um)  and 
are  shown  in  Figure  IV-20.  Consider  first  the  upper  dashed  cur\e. 

It  is  the  spectral  flux  density  (watts/cm  of  the  residual  background 

O 

signal  if  the  temperature,  of  the  particle  surroundings  were  300  K 

o 

and  the  temperature,  T^^,  of  the  chopper  blade  were  301  K.  The  next 

O O 

lower  dashed  cur\'c  is  for  Tj^  = 300  K and  T^  = 300.2  K,  a temperature 

O O 

difference  of  0.  Z K,  The  lowest  dashed  curve  Is  for  T,,  = 80  K and 

® r I 1 

Tp  = 85  K.  The  solid  curves  are  the  spectral  flux  density  (watts/cm  ) 
collected  by  the  particle  emissivity  apparatus  from  a black  Z^m  radius 

O 

particle  for  the  indicated  temperatures  T . For  the  case  T = 1000  K 
and  Tjj  - Tj^  = 0.  2 K the  particle  signal  is  comfortably  above  the  back- 
ground signal.  However,  for  a particle  whose  temperature,  T^  = 500  K, 
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the  background  signal  Is  significantly  higher  than  the  particle  signal  up 
to  ISOOcm"^  (•'-  6.  6|j.m).  As  long  as  T is  at  least  a few  hundred  K 
greater  than  Tg  and  is  close  to  Tg,  the  particle  emission  spectrum 
can  be  obtained  by  differencing  spectra  obtained  with  and  without  the 

O 

heated  particle  present.  For  particle  temperatures  much  below  1000  K 
this  means  differencing  two  numbers  of  comparable  magnitude  and  thus 
accuracy  is  lost.  In  fact,  as  T approaches  Tg  this  difference  goes  to 
zero  and  no  data  can  be  obtained.  The  several  curves  shown  in 
Figure  IV-20  suggest  the  specific  values  of  and  Tg  - which  must 
be  achieved  for  the  experiment  to  produce  meaningful  data. 


b.  Particle  Temperature  and  Emissivity  Determination 

The  determination  of  particle  emissi\  ity  from  the  observed 
radiation  spectra  has  been  recognized  as  a critical  analytical  step  since 
the  inception  of  the  POPM  program  and  has  therefore  been  the  subject 
of  continued  analysis.  The  determination  of  the  emissivity  implies  a 
knowledge  of  the  particle  temperature,  which  can  be  obtained  most 
simply  from  an  analysis  of  the  infrared  radiation  emitted  by  the 
particle. 

The  statement  of  work  for  the  POPM  program  dated  14  June  1974 
contains  a proposed  analysis  scheme  for  obtaining  temperatures  and 
emissivities  from  considerations  of  the  spectral  \ariation  of  the  ratio 
of  spectra  of  the  same  particle  obtained  at  two  different  temperatures. 
Subsequent  analysis  has  shown  that  this  procedure  was  very  dependent 
on  data  at  500cm  ^ and  below,  and  that  small  errors  in  spectral 
intensity  values  lead  to  large  uncertainties  in  the  temperature. 

A new  metliod  to  determine  particle  temperature  and  emissi%  ity 
based  on  auxiliary  observations  of  a blackbody  at  a known  high  temper- 
atures has  therefore  been  de\  ised.  This  subsection  presents  the 
details  of  this  data  analysis  method. 
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(1)  Particle  Emlsslvlty  Determination 

Neglecting  the  background  contribution  (for  which  we 
shall  be  able  to  correct),  the  intensity  of  radiation  as  seen  by  a spec- 
troscopic device  (in  our  case  a k.  ichelson  interferometer)  can  be 
written. 


Jp(v,  Tp)  dv 


Tie  (v,T  ) Cjr^v^f(v)  dv 


hcv 

kT 

. P . 1 


where 


e (v,  T)  = emissivity  of  the  particle  which  in 
^ general  may  be  a function  of  v,  and  T, 

Cj  = first  radiation  constant, 

Tp  = radius  of  the  particle, 

f(v)  = function  of  v accounting  for  instrumental 

transmission,  detector  efficiency,  and 
atmospheric  absorption, 

h = Planck's  constant, 

c = velocity  of  light, 

k = Boltzmann's  constant, 

V = wavenumber  of  radiation  (=  X"^),  and 

T = particle  temperature. 


For  a blackbody  radiation  source  a similar  equation  may  be 
written  except  by  definition  T)  = 1,  and  in  practice  manufac- 

turers claim  to  achieve  6gg(v,  T)  > 0.  995.  The  blackbody  radiance 
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nC^r|v^f(v)  dv 


where  r^  is  the  radius  of  the  circular  aperture  of  the  blackbody  radi- 
ating to  the  outside  world. 

Dividing  equation  IV-4  by  equation  IV- 5 we  obtain 


e (v,  T ) 
P P 


hcv 

3 ^ - 1 

hcv 

kT 

. P - 1 


J and  J will  be  measured  experimentally,  r„  and  are  know\. 
p XD  i5 

The  radius  of  the  particle  r will  be  determined  experimentally  by 
measuring  the  scattered  light  as  a function  of  scattering  angle. 
Solving  for  e (v,  T ), 

Jr 


e (v,  T ) 
P P 


"b^p 


Note  that  by  dividing  a particle  emission  spectrum  by  a blackbody 
emission  spectrum  we  obtain  an  expression  for  Sp('^«  T^)  directly, 
providing  Tp  can  be  obtained. 

The  accuracy  of  the  determination  of  Sp  depends  critically  on 
the  determination  of  the  particle  temperature.  Let  the  error  in  Tp 
be  ‘^Tp,  then 
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he  V 


and  using  equation  IV-7 


he  V 


Thus  the  fraetional  error  in  e equals  the  fraetional  error  in  T 

P X X P 

multiplied  by  a funetion  of  the  form  xe  /(e  - 1),  where  x = hev/kT  . 

P 

The  minus  sign  indicates  the  errors  run  counter  to  each  other.  If 

T is  overestimated,  e will  be  underestimated  and  vice  versa.  Note 
P P 

that  the  relationship  is  a function  of  x/T  alone.  In  the  particle  emis- 

-1  P 

sivity  experiment  ~ 400cm  (X  - 25|a.m)  and  T^lmax)  ~2Z00°K. 

Thereforex  ~0.26.  Aplotofxversusxe^/(e^-l)outtox=7 

is  shown  in  Figure  IV-21.  As  x approaches  zero  xe^/(e^  - 1)  goes  to  1. 

Thus,  for  a given  percentage  error  in  T the  error  in  e for  small 

P P 

v/T  will  be  of  the  same  order  as  fiT  /T  . From  the  graph  Table  IV-4 
P P P 

was  constructed  for  a one  percent  error  in  T at  various  temperatures 

P 

and  radiation  frequencies.  The  above  analysis  should  be  used  with 

caution  in  that  the  error  propagated  into  Sp  by  an  error  in  T may  not 

always  be  the  most  important.  For  instance,  even  a "black"  500  °K 

particle  will  not  radiate  much  power  in  the  5000cm"  ^ (\  = 2um)  region, 

and  therefore  the  determination  of  Sp  may  be  governed  more  by  how 

well  Jp  is  determined  rather  than  T^.  However,  when  the  S/N  in  J 

is  sufficiently  high,  the  error  in  the  determination  of  Gp  will  be 

largely  governed  by  the  error  in  T . Therefore,  even  if  J has  a 

P P 

small  uncertainty,  fractional  errors  in  T will  result  in  fractional 

P 

errors  from  two  to  ten  times  larger  in  e . 

P 


Table  IV-4.  Error  In  e for  a one  percent  error  In  T 


500cm 


1 500cm 


2500cm 


5000cm 


(2)  Determination  of  T 

P 

The  reliable  determination  of  T is  central  to  the 

P 

successful  determination  of  Sp(v,  T^).  If  Sp(v,  T ) is  reasonably  well 
behaved,  (the  meaning  of  "reasonably  well  behaved"  will  become 
apparent  further  in  this  section)  in  one  or  more  fairly  broad  wave- 
number  regions  for  which  v/T^  s 3.  5(cmK)'^,  then  T can  be  deter- 
mined. 

Recalling  equation  IV-6,  and  putting  all  the  knowns  on  the  left 
side  of  the  equation,  we  can  define  a K(v) 


‘^p^B  / hcv 


le 


e (v,T  ) 
P P 


he  V 
kT 

; P - 1 


All  the  quantities  on  the  left  are  known,  therefore  Kfv)  is  known.  Now- 

in  the  high  wavenumber  end  of  the  spectrum,  where  v/T  s 3.5(cmK)'^ 
^ fhev"!  ^ , P 

so  that  exp  » 1,  we  have 

- P- 

-hev/kT 

K(v)  = e (v,T  )e  P,  or 

P P 

(I 

InK(v)  = Ine  (v,T)-{j^ 

^ P 


If  e is  a constant  over  some  wavenumber  inter\  al  for  which 

P - 1 

'-/Tp  3.  5(cmK)  , then  a plot  of  InK(v)  versus  v will  be  a straight 

line  whose  slope  is  - hc/kT  , thus  T will  be  determined.  Obviously, 

P P 

the  wider  the  interval(s)  the  more  accurately  Tp  will  be  dete  rniined. 

IZxcept  in  the  region  of  sharp  absorption  bands  (and  therefore 
probably  emission  bands)  the  emissivity  would  be  expected  to  be  a 
slowly  varying  function  of  v.  For  example,  let  us  assume  over  some 
wavenumber  intersal 
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i 


I 


ep(v,  Tp)  = a(l  - bv) 


such  that  b Is  small  enough  so  that 


e (v,  T ) = a(l  - bv)  ae 
P P 


Ine  = Ina  - bv 
P 


and  equation  IV-8  would  become 


InK(v)  = Ina  - bv  - 


, hcv  / 1 bk 


Obviously,  b must  be  very  much  less  than  T if  T is  to  be 

P P o 

determined  accurately.  For  example,  If  T lies  In  the  range  800  K 

o ^ - 1 

to  1500  K and  Sp  ~ 0.  3 (average  in  a 1000cm  range  somewhere 
between  4000  and  8000cm  ^),  variation  of  £p  of  three  to  six  percent 
as  given  by  equation  IV-7,  would  still  allow  Tp  to  be  determined  to 
about  one  percent.  Particle  emissivities  have  not  been  measured, 
but  bulk  emissivities  of  the  materials  of  interest  show  approximately 
this  variation. 


Consider  now  the  emission  spectrum  of  a particle  measured  at 
two  temperatures  Tp^  and  Tp^.  Taking  the  ratio  of  these  spectra  we 
obtain 


J (v,T  ,) 

p pi 


P pi 


hcv 
kT  , 

> - 1 
hcv 
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and  for  v/T  , and  v/T  ^ s 3.  5(cmK)  ^ we  have 
pi  p-i 


hc\ 


■pi  p2j 


and  therefore 


In 


Jp("' 


= In 


e (v,  T ,) 

P .P^.. 

e (v,  T , ) 
p'  pi' 


hcv 


•pi  "p2y 


(lV-11) 


If  the  ratio  of  the  emlsslvitles  Is  a constant  then  the  above 
expression  Is  the  equation  for  a straight  line  whose  slope  Is 

he  / 1 1 

k ( T " T , 

V pi  P2 

(3)  Particle  Radiance  Data  Reduction 

The  preceding  analysis  suggests  a procedure  for 
obtaining  particle  temperatures  from  the  particle  emission  spectra. 
We  will  obtain  spectra  at  a number  of  unknown  temperatures  and 
blackbody  emission  spectra  at  a number  of  known  temperatures.  It 
appears  that  four  temperatures  will  be  adequate;  ~ 2000,  1600,  1200 
and  800  °K.  The  slgnal-to-nolse  ratios  for  these  spectra  are  antici- 
pated to  be  much  greater  than  for  the  particle  emission  spectra.  Rho 
(p)  curves  of  the  type 


PpB 


(particle) 

(blackbody) 
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will  be  generated  and  Inc  _ versus  v will  be  examined  in  the  2000  to 

!]  P® 

4000cm  interval  for  straight  line  behavior.  In  a similar  fashion 
c curves  of  the  type 

J(particle, 

^PP  ~ J(particle,  T^) 

will  be  generated  and  examined  for  straight  line  behavior  ‘"“pp 
versus  V.  In  order  to  simplify  the  notation,  define 

X.  = reciprocal  of  the  temperature  of  a blackbody  at 
temperature  T.^,  and 

Y = reciprocal  of  the  temperature  of  a particle  at 

j 

temperature  T.. 

Each  physically  acceptable  straight  line  plot  based  on  two 
spectra  results  in  a slope,  proportional  to  the  difference  between 

the  reciprocals  of  the  temperatures  associated  with  the  two  spectra 
for  which  the  logarithmic  ratio  has  been  plotted.  The  result  of  all 
the  possible  spectral  ratios  is  a system  of  equations  of  the  form: 


X.  - Y.  = M.. 
i I tj 


Y,  - Y = M, 
k e ke 

The  X's  and  M's  are  known,  and  the  Y's  are  to  be  determined.  Si 
more  than  one  segment  of  the  spectra  may  be  employed,  possibly 
yielding  multiple  \alucs  of  M for  the  same  spectra,  the  Y's  will 


The  noise  level  in  the  experimental  data  is  determinable  (see 
Section  IV-B.4.a),  Criteria  can  tlierefore  l)c  set  up  to  determine  if 
In'  is  a straight  line  within  the  known  signal-to-noisc. 

'pB 
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probably  be  over  determined.  This  can  be  combined  with  additional 
a priori  knowledge  to  Improv^e  the  accuracy  of  the  Y's  (and  thus  tem- 
peratures determined).  The  ordering  of  the  particle  temperatures 
will  be  known  since  the  highest  laser  power  should  correspond  to  the 
highest  temperature,  the  second  highest  laser  power  to  the  second 
highest  temperature,  and  so  on.  Equations  of  the  type  of  IV-12  deter- 
mine the  Y's  temperature.  Equations  of  the  type  of  IV-13  can  be  used, 
along  with  the  ordering  information  from  laser  powers  to  determine  if 
the  correct  slope  value  has  been  used  in  the  type  IV-12  equations. 


! 


Once  Tp  for  a particle  emission  spectrum  is  determined,  the 
emissivity  can  then  be  obtained  from  the  ratio  of  the  particle  spec- 
trum to  a suitable  blackbody  spectrum  times  the  calculable  expres- 
sions shown  in  equation  IV- 7. 


Samples  for  Particle  Emissivity  Studies 


The  materials  and  size  ranges  for  particles  of  interest  for 
emissivity  measurements  in  the  present  program  are:  Carbon, 
having  a diameter  from  0.  1 to  l.Ogm;  magnesia  (MgO),  zirconia 
(ZrO^)>  and  alumina  (Al^O^)  having  diameters  from  1 to  25-m.  The 
study  was  restricted  to  spherical  particles  to  facilitate  the  compari- 
son of  measured  emissivity  with  computed  values. 


a.  Materials 

(1)  Carbon 

Pure  solid  carbon  exists  in  three  stable  forms. 
The  two  crystalline  forms  are  graphite  and  diamond  and  the  third  is 
ouasi-g raphitic  and  known  as  carbon  black.  Our  main  concern  here 
is  with  carbon  black  since  that  is  the  form  of  material  loaded  into 
rocket  motors  and  most  probably  the  form  in  which  it  is  present  in 
rocket  exhausts. 


-208- 


Carbon  black  In  the  size  range  of  interest  for  the  particle  emis- 
sivity  experiment  is  produced  by  the  thermal  decomposition  of  natural 
gas  in  the  absence  of  air  or  flame  and  is  called  thermal  black.  Carbon 
blacks  display  two-dimensional  crystallinity.  They  consist  of  well 
developed  graphite  platelets  stacked  roughly  parallel  to  one  another  but 
random  in  orientation  about  the  normal  to  the  layers.  When  carbon 
black  is  heated,  the  growth  of  the  quasi-graphitic  parallel  layer  groups 

o 

begins  at  about  1000  C.  In  the  larger  carbon  black  powders  such  as 
thermal  black,  three-dimensional  graphitic  structure  is  detected  after 

O 

heating  to  3000  C.  In  the  finer  carbon  blacks,  a three-dimensional 
graphitic  structure  apparently  never  evolves.  Thermal  blacks  have 
negligible  oxygen  content,  and  hydrogen  present  is  from  the  hydro- 
carbon raw  material  and  is  distributed  throughout  the  particle.  The 
finer  carbon  blacks  can  have  appreciable  oxygen  present  on  the  sur- 
face of  each  particle.  The  oxygen  and  hydrogen  combine  with  surface 
carbon  atoms  and  form  a variety  of  surface  complexes.  On  heating 
carbon  black  progressively  from  600  to  1500°C,  decomposition  pro- 
ducts of  these  surface  compounds  are  evolved  as  CO^,  CO  and  H^. 

In  air,  serious  oxidation  of  carbon  probably  begins  around 

O 

500  C.  This  probably  does  not  effect  the  surface  chemistry,  however, 
since  all  the  products  of  oxidation  are  gaseous.  Oxidation  will  cause 
a small  particle  to  eventually  be  converted  entirely  to  the  gas  phase. 

(2)  Magnesia 

The  chemical  properties  of  magnesium  oxide  vary 
considerably  with  the  purity  of  the  niaterial.  MgO  prepared  in  the 
temperature  range  from  400  to  900  C rehydrates  readily  e\en  in  cold 
water.  Fused  magnesium  oxides  are  virtually  insoluble  in  water 
except  in  extremely  fine  powders.  Because  of  the  rapid  hydration 
brought  about  by  the  presence  of  impurities,  it  is  essentia]  to  oistain 
as  pun-  MgO  as  possible. 


f 
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MgO  crystallizes  from  the  liquid  in  the  face  centered  cubic 
system.  As  with  Al^O,  there  is  only  one  stable  crystal  form  of  pure 
MgO  in  the  range  from  room  temperature  to  the  melting  point  at 
ZBOO'C. 

(3)  Zirconia 

ZrO^  exists  stably  at  ordinary  temperatures  as 
monoclinic  crystals.  Transformation  between  monoclinic  and  tetra- 
gonal symmetry  occurs  at  about  1200 °C.  This  transformation  is 
accompanied  by  a large  change  in  volume,  about  7.  5 percent.  This 
gives  rise  to  cracks  and  failures  in  crystals  and  in  ceramic  bodies 
cycled  through  this  temperature.  This  can  be  circumvented  by 
inducing  the  zirconia  to  crystallize  in  the  cubic  system,  through  the 
influence  of  certain  dissolved  oxides. 

All  natural  zirconium  ores  contain  hafnium.  Concentrations 
vary  from  about  two  percent  to  somewhat  more.  The  process  of 
removing  this  hafnium  is  expensive,  and  the  zirconium  dioxide  powder 
used  to  load  rocket  motors  most  probably  contains  hafnium.  Impure 
zirconia  costs  about  $2.  15/lb,  while  zirconia  of  comparable  size  and 
purity,  but  with  only  25ppm  of  hafnium  costs  about  $500/lb.  There 
does  not  seem  to  be  any  compelling  need  to  use  this  costly  grade;  so 
the  ZrO^  materials  studied  in  this  experiment  are  designated  "pure”, 
when  in  fact  they  contain  two  percent  hafnium. 

(4)  Alumina 

a-Al^O^  crystallizes  from  the  liquid  in  the  rhom- 
bohedral  system.  This  is  also  the  only  stable  form  of  pure  anhydrous 
Al^Oj  at  all  temperatures  from  room  temperature  to  the  melting  point 
at  2040  C,  Because  of  this,  there  is  no  difficulty  in  insuring  that  all 
forms  of  pure  solid  Al^O^  will  have  the  same  crystal  structure  in 
these  experiments.  It  is  produced  industrially  by  solidification  of 
molten  alumina,  from  a sintering  process  in  the  temperature  range 
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1 rotn  1 300  to  1700  C,  and  also  by  calcination  at  a sufficiently  high 
temperature  (1100  to  1300  C)  of  the  other  forms  of  alumina.  There 
arc  a great  niany  stable  forms  of  aluminum  oxide  when  mixed  with 
hydroxyl  groups  to  produce  the  so-called  hydrated  aluminas.  It  is 
necessary,  therefore,  to  obtain  as  pure  a form  of  A1,0_  as  possible 
since  a number  of  impurities  can  cause  rapid  hydration  of  a-Al^O^. 

b.  Spheroidal  Sample  Preparation 

(1)  Carbon 

Carbon  powders  are  available  as  amorphous  carbon, 
graphite  and  dianiond.  Diamond  is  not  available  in  spherical  form  nor 
is  it  possible  to  produce  it  in  that  form.  This  results  from  the  funda- 
mental crystal  structure  of  this  material  which  imposes  flat  surfaces 
on  the  particle.  The  primary  interest  of  this  program,  in  any  case, 
i is  in  amorplious  carbon,  that  being  the  form  of  the  material  most 

I probably  present  in  rocket  exhaust. 

The  powders  of  spherical  particles  that  we  ha\e  obtained  are  of 
very  high  purity  and  meet  the  size  requirements.  The  Kreka  Corpor- 
ation supplied  us  with  amorphous  carbon  particles  that  are  almost 
perfect  spheres  in  the  1 . 5 to  Sum  range  and  are  99.  9 percent  pure 
(Figure  IV-22).  For  the  0.3  to  0.  Sum  particle  size  range,  the  Huber 
and  Vanderbilt  Corporations  supplied  us  with  powders  that  arc  mostly 
spherical  and  are  99  percent  pure. 

(2)  Metal  Oxides 

Abundant  sources  of  these  materials  in  powder 
forms  are  available  in  a wide  assortment  of  sizes.  Howe\er,  no 
suppliers  of  any  of  these  materials  can  claim  capability  of  pro\  iding 
! spherical  particles.  Conversations  with  manufacturers  of  plasma 

arc  heaters,  on  the  other  hand,  revealed  that  it  is  possible  to  form 
spherical  particles  by  melting  the  powders  in  a plasma  spray. 


Plasmadyne  Corporation  of  Anaheim,  California  has  success- 
fully generated  spheres  of  AI2O2  and  ZrO^  in  the  size  range  from  20 
to  80um.  Photomicrographic  studies  of  these  alumina  and  zirconia 
powders  show  that  while  spheres  do  indeed  exist,  the  yield  of  spheres 
is  only  about  ten  percent.  Apparently,  this  process  functions  in  such 
a way  as  to  con\ert  particles  to  spheres  in  only  over  a fairly  narrow 
range  of  particle  masses.  Larger  particles  do  not  melt  completely 
and  smaller  particles  evaporate  and  subsequently  recondense  as 
many  very  small  particles.  Plasmadyne  has  not  been  successful  in 
converting  MgO  to  spheres.  Tafa  Industries  of  Bow,  New  Hampshire 
claim  to  be  able  to  produce  spheres  in  the  size  range  of  2 to  20um 
from  any  one  of  the  three  oxides  of  interest.  Therefore,  it  would 
seem  that  if  these  companies  were  supplied  with  the  proper  sized 
powders  that  spherical  particles  could  be  obtained. 

The  Aerospace  Corporation  is  capable  of  supplying  these  com- 
panies with  a narrow  range  of  particle  sizes,  as  we  are  equipped  with 
a Donaldson  particle  classifier  that  separates  any  range  of  particles 
desired  from  powders  with  a large  size  distribution. 

Several  sources  for  these  oxides  were  contacted,  photomicro- 
graphic studies  of  many  samples  were  made  and  the  powders  that 
could  meet  the  program  requirements,  after  classification  and 
s phe roidization  were  chosen. 

(a)  Magnesia 

We  have  procured  three  powders  that  together 
meet  all  of  the  requirements  of  the  experiment  after  classification 
and  spheroidization:  Fisher  Chemical  M51 -75617,  Kaiser  Chemical 
325  mesh,  Ventron  325  mesh.  No  attempt,  as  yet,  has  been  made 
to  classify  or  spherisize  these  particles. 
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(b)  Zirconia 

Two  powders  have  been  located  which,  after 
classification  and  spheroidizatlon  will  meet  the  requirements:  Fisher 

Chemical;  2-83,  ML/B,  2x120/8404.  No  attempt,  as  yet,  has  been 
made  to  classify  or  spherisize  these  particles. 

(c)  Alumina 

The  size  range  from  7-9^im  of  the  Fisher 
alumina  A-591  was  separated  using  the  classifier  and  sent  to  Tafa 
Corporation  for  spheroidizatlon.  This  attempt  failed;  the  powder 
agglomerated  and  could  not  be  plasma  sprayed.  However,  by  using 
their  own  Al^O^  powders  as  a starting  material,  Tafa  Corporation 
was  able  to  supply  Aerospace  with  spherical  particles,  "lonarc” 
Al^O^-3ZO  and  400.  These  particles  have  a size  distribution  from 
4-20um  and  15-30^im,  respectively,  and  are  shown  in  Figure  IV-23. 


C.  SUMMARY,  DISCUSSION  AND  RECOMMENDATIONS 

Particles  of  all  four  materials  of  interest  in  this  study  have  been 
levitated  and  maintained  in  place  while  being  irradiated  with  intense 
1 , 06u,m  Nd: YAG  laser  radiation. 

A method  of  analyzing  the  spectral  emission  data  has  been 
derived  and  is  presented  in  this  report.  The  emissivity  depends  on 
finding  the  particle  temperature.  A procedure  for  doing  this  is  pre- 
sented which  depends  on  finding  a spectral  region,  v>  2000cm  \ in 
which  the  emissivity  does  not  vary  rapidly  with  frequency  over  a 
sufficient  frequency  interval  that  a slope  can  be  determined.  In  view 
of  the  computed  particle  emissivities  shown  in  Section  III,  this  condi- 
tion may  not  be  met  for  some  of  the  materials  of  interest. 

Samples  of  all  the  materials  are  available  in  the  size  range 
required  for  the  particle  experiments.  Spherical  particle  samples  of 
Al^O^  and  carbon  are  available. 
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Due  to  the  problems  encountered  in  simultaneously  levitating 
and  irradiating  particles,  no  infrared  radiation  has  been  observed 
from  the  hot  particles  within  the  period  covered  by  this  report. 

We  believe  that  significant  evidence  exists  that  the  experiment 
is  possible  and  recommend  that  it  be  continued.  Specifically,  we 
recommend  the  following  approach: 

(1)  Use  an  infrared  detector  with  the  required  collection 
optics  to  observe  the  infrared  radiation  from  the  hot 
particle.  This  will  provide  a demonstration  of  the  feasi- 
bility of  the  overall  experiment. 

(2)  Measure  the  particle  emission  spectrum  with  a Michelson 
interferometer  spectrometer.  The  optics  should  be  con- 
figured so  that  calibration  spectra  from  a high  temperature 
blackbody  can  also  be  obtained.  The  emissity  spectrum 
will  then  be  deduced  by  the  methods  outlined  in  this  section, 
or  an  improvement  of  these  procedures. 
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SECTION  V 


SUMMARY  AND  RECOMMENDATIONS 


Tlu'  o\aluation  of  the  reliability  of  using  bulk  refractive  index 
data  to  calc  ulate  infrared  particle  emissivity  has  been  the  basic 
motivation  cor  the  POPM  program.  In  this  report  we  have  shown 
that  particle;  emissivity  calculations  from  bulk  data  may  indeed  be 
\ery  unreliable.  This  emphasizes  the  need  for  continued  efforts  to 
directly  determine  infrared  particle  emissivity.  To  this  end  an 
apparatus  has  been  assembled  which  levitates  and  simultaneously 
heats  single  particles.  Analytical  means  have  been  derived  to 
determine  encissic  itles  from  experimental  spectral  radiance  meas- 
urements. The  significant  achievements  in  each  of  the  areas  covered 
in  this  report  and  leading  to  the  general  conclusions  stated  above  are 
summarized  in  this  section.  For  a more  detailed  discussion  consult 
the  appropriate  section  of  this  report. 


A.  lU’ LK  EXPERIMENT 

A de\ice  was  designed  and  constructed  for  the  measurement  of 
the  spectral  reflectance  of  polarized  infrared  radiation  by  plane  sur- 
faced samples  at  room  temperature,  and  higher,  fc^r  a range  of  inci- 
dence angles.  The  originally  planned  inversion  of  the  Fresnel  equations 
for  various  angles  to  obtain  the  complex  refractive  index  was  found  to 
provide  ambiguous  results.  The  refractive  index  data  were  therefore 
obtained  by  either  classical  dispersion  analysis,  or  Kramers-Kronig 
analysis,  of  the  reflectivity  data  obtained  at  the  fixed  incidence  angle 
of  45  degrees. 

Reflectivity  data  between  ~ 4.  5 and  25gm  have  been  obtained  and 
analyzed  to  determine  n and  k for  Al^O^,  MgO  and  ZrO^  at  room  tem- 
perature and  a temperature  '^600°K.  Similar  measurements  for 

O 

carbon  have  been  made  only  at  300  K. 
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In  general,  the  present  refractive  index  data  fall  within,  or  near 
the  range  of  values  given  in  the  literature,  which  are  mostly  room 
temperature  values.  The  values  at  higher  temperatures  behave  quali- 
tatively as  expected,  the  reststrahl  lattice  frequencies  moving  slightly 
to  lower  frequencies,  with  an  increase  in  bandwidth. 

As  error  analysis  of  the  carbon  data  has  been  made  indicating 
errors  of  five  to  ten  percent  in  n and  k for  most  of  the  refractive  index 
data,  and  higher  uncertainties  for  other  values,  notably  for  k < 0.  05. 
Direct  reflectivity  measurements  are  not  well  suited  to  determine 
small  values  of  k.  As  discussed  previously,  transmissivity  measure- 
ments or  total  attenuated  reflectivity  measurements  are  capable  of 

_4 

measuring  k down  to  10  , or  less. 

The  temperatures  for  which  data  are  presented  here  are  below 
those  thought  typical  of  particles  in  applications.  The  measurements 
were  limited  by  sample  heater  failure,  which  we  do  not  belie\e  to  be 
a major  limitation. 

The  continued  bulk  measurement  program  should  incorporate 
the  following  improvements  to  the  basically  operational  bulk  experi- 
ment. 

(1)  Improve  the  thermal  contact  between  the  heater  and 
sample  so  that  sample  temperatures  in  excess  of  1000°K 
may  be  achieved. 

(2)  Improve  the  chopping  system  in  the  experiment  to  provide 
more  efficient  chopping.  This  should  improNC  the  signal - 
to-noise  ratio  by  at  least  a factor  of  three  to  five  and 
thus  lower  the  errors  in  reflectivity  nieasurcnients. 

(3)  Determine  theoretically  if  45  degrees  is  the  optimum 
incidence  angle  for  refractive  index  determination. 

(4)  Devise  a technique  complementary  to  the  present  reflect- 
ance technique  to  determine  accurately  low  values  of  k, 
especially  for  heated  samples. 


-220- 


B. 


PARTICLE  EMISSIVITY  COMPUTATIONS 


The  emissLvLtles  expected  for  spherical  particles  of  a variety 
of  radii  have  been  calculated  by  Mie  theory  from  complex  refractive 
indices  determined  in  the  bulk  refractive  index  measurements  pro- 
gram. These  calculations  indicate  the  emissivity  of  small  particles 
depends  critically  on  parameters  which  are  very  likely  to  change 
markedly  between  bulk  and  particle  conditions. 


The  calculations  show  emissivities  greater  than  unity  which  are 
caused  by  expected  resonance  effects  between  vibration  modes  of  the 
particle  and  the  infrared  radiation. 

The  emissivities  are  critically  dependent  on  size. 

The  emissivities  are  strong  functions  of  frequency  in  the  spectral 
range  corresponding  to  the  reststrahl  bands  in  the  bulk  material.  In 
this  same  spectral  region  the  interaction  between  particle  size  and 
refractive  index  is  strong,  complicated  and  difficult  to  predict  without 
a full  Mie  theory  analysis  for  particles  in  the  resonant  frequency  size 
range.  The  reader  is  reminded  of  the  Al^O^  particle  emissivity 
spectra  presented  in  Section  III  as  an  example. 

In  the  spectral  region  of  wavelength  short  with  respect  to  parti- 
cle size  and  where  low  values  of  k occur,  the  emissivity  depends 
directly  on  the  value  of  k.  See  Figures  III-27  and  111-28  for  example. 


The  observed  emissivity  of  particles  large  compared  with  the 
wavelength  of  radiation  and  away  from  the  reststrahl  regions  is 
expected  to  be  different  than  the  computed  emissivity  because  of  the 
dependence  of  emissivity  on  k and  the  liklihood  that  \-alues  of  k in 
particles  will  differ  significantly  from  bulk  values  due  to  imperfec- 
tions in  the  particles.  This  result  of  the  theoretical  calculations 
emphasizes  the  importance  of  observing  experimentally  the  emission 
from  single  particles. 

The  theoretical  analysis  effort  should  be  continued  to  determine 
the  effects  of  particle  non- sphericity  and  surface  roughness  on  the 
emissivity  of  particles. 


-221- 


C.  PARTICLE  EMISSIVITY  EXPERIMENT 

An  experimental  apparatus  has  been  designed  and  constructed 
which  has  demonstrated  the  capability  to  levitate  and  simultaneously 
irradiate  with  Nd:YAG  laser  radiation  single  particles  of  Al^O^.  C, 

MgO  and  ZrO^*  Analysis  indicates  the  observed  laser  intensities 
should  be  sufficient  to  heat  the  particles  to  temperatures  approaching 
the  melting  point  of  Al^O^. 

Analytical  procedures  have  been  devised  to  determine  the 
emissivity  of  the  particles  from  observed  spectral  measurements  of 
the  intensity  of  the  emitted  radiation,  providing  a spectral  region  of 
nearly  uniform  particle  emissivity  can  be  found  in  the  2000-8000cm"  ^ 
region. 

The  difficult  problem  of  simultaneously  levitating  and  irradiating 
a particle  was  solved  shortly  before  the  close  of  the  time  period 
covered  by  this  report.  Therefore,  no  thermal  energy  from  a single 
particle  has  been  observed. 

The  continued  single  particle  emissivity  experiment  should 
incorporate  the  following  steps: 

(1)  Detect  thermal  radiation  from  a laser  heated  levitated 
particle.  This  will  provide  an  empirical  demonstration 
of  the  feasibility  of  the  experiment. 

(2)  Incorporate  the  fast  scanning  Michelson  interferometer 
and  high  temperature  calibration  source  into  the  experi- 
mental apparatus  to  measure  the  spectral  dependence  of 
the  emitted  radiation. 

(3)  Determine  from  the  emission  measurements  the  emissivity 
of  small  spherical  particles  and  compare  with  the  theoreti- 
cally computed  emissivity  values. 

(4)  Extend  the  measurements  on  single  particles  to  non-spherical 
shapes  and  other  materials. 
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GENERAL  CONCLUSIONS  AND  RECOMMENDATIONS 


i 


D. 

The  computed  emisslvity  of  small  particles  can  be  expected  to 
differ  widely  from  the  observed  emisslvity,  based  on  computations 
presented  in  Section  III  and  summarized  above.  This  provides  a 
partial  answer  to  the  fundamental  question  to  which  the  POPM  pro- 
gram has  been  addressed.  This  analytically  derived  answer  should 
be  verified  by  empirical  single  particle  observations. 

We  recommend  that  in  the  near  future  the  highest  priority  in 
the  POPM  program  should  be  given  to  the  single  particle  emisslvity 
experiment.  A lower  priority  should  be  given  to  the  theoretical 
analysis  of  non-spherical  shapes.  At  this  time  the  bulk  refractive 
index  is  also  of  low  priority,  but  in  the  more  distant  future  bulk  data 
will  be  required  to  analyze  the  then  available  single  particle  data. 
Therefore  the  improvements  in  this  bulk  experiment  outlined  above 
should  be  implemented  on  the  basis  of  non-interference  with  the 
particle  experiment. 


I 
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THE  IVAN  A.  GETTING  LABORATORIES 

Th«  Laboratory  Oparationa  of  Tha  Aaroapaca  Corporation  ia  conducting 
axparimcntal  and  thaoratical  invaatigationa  nacaaaary  for  tha  avaluation  and 
application  of  acientific  advancea  to  new  military  cencepta  and  ayatema.  Ver- 
aatillty  and  flasibility  hava  baan  developed  to  a high  degree  by  the  laboratory 
paraonnal  in  dealing  with  the  many  problama  encountered  in  the  nation'a  rapidly 
developing  apace  and  miaaile  ayatema.  Expertiae  in  the  lateat  acientific  devel- 
opmenta  ia  vital  to  the  accompliahment  of  taaka  related  to  theae  problema.  The 
laboratoriaa  that  contribute  to  thia  reaearch  are: 

Aerophvaica  Laboratory:  Launch  and  reentry  aerodynamica,  heat  trana* 
far,  reant^  pkyaica,  chemical  kinetica,  atructural  mechanica,  flight  dyaamica, 
atmoapheric  pollution,  and  high>power  gaa  laaera. 

Che^atry  and  Phyaica  Laboratory:  Atmoapheric  reactiona  and  atmoa- 
pheric~^tlca,  chemical  reacttona  in  polluted  atmoapherea,  chemical  reactiona 
of  excited  apeciea  in  racket  plumea,  chemical  thermodynamica,  plaama  and 
laaer-induced  reactiona,  laaer  chemiatry,  propulaion  chemiatry,  apaca  vacuum 
and  radiation  effecta  on  materiala,  lubrication  and  aurface  phenomena,  photo- 
aenaitive  materiala  and  aenaora,  high  precialon  laaer  ranging,  and  tha  appli- 
cation of  phyaica  and  chemiatry  to  problema  of  law  enforcement  and  biomedicine. 

Electronica  Reaearch  Laboratory:  Electromagnetic  theory,  devicea,  and 
propagation  phenomena,  including  plaama  electromagnetica:  quantum  electronica, 
laaara,  and  electro-optica;  communication  aclencea,  applied  elactronica,  aemi- 
condueting,  auperconducting,  and  cryatal  device  phyaica,  optical  and  acouatlcal 
imaging;  atmoapheric  pollution;  millimeter  wave  and  far-infrared  technology, 

Materiala  Sciencea  La^ratorv:  Development  of  new  materiala;  metal 
matrix  compoaitea  and  new  /orma  or  carbon;  teat  and  evaluation  of  graphite 
and  ceramica  in  reentry;  apacecraft  materiala  and  electronic  compouenta  in 
nuclear  weapona  environment;  application  of  fracture  mechanica  to  atreaa  cor- 
roaion  and  fatigue-induced  fracturea  in  atructural  metala. 

Space  Sciencea  Laboratory:  Atmoapheric  and  ionoapheric  phyaica,  radia- 
tion train  tha  atmoaphere,  daneity  and  compoaitioa  of  the  atmoaidiare,  aurorae 
and  oirglow;  magnataapheric  phyaica,  coai^c  raya,  generation  and  propagation 
of  plaama  wavea  in  the  magnetoaphere;  aolar  phyaica,  atudlea  of  aolar  magnetic 
flelda;  apace  aatronomy,  x-ray  aatronomy;  the  effecta  of  nuclear  exploaiona, 
magnetic  atorma,  and  aolar  activity  on  the  earth'a  atmoaphere,  ionoepherc,  and 
magnetoaphare;  the  effecta  of  optlcaL  electromagnetic,  and  particulata  radia- 
Hona  la  apace  on  apace  ayatema. 
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